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GRAVIMETRIC GEODESY AND SEA SURFACE TOPOGRAPHY STUDIES 
BY MEANS OF SATELLITE-TO-SATELLITE TRACKING 
AND SATELLITE ALTIMETRY 
Joseph W. Siry 

ABSTRACT 

The present knowledge of the earth's gravitational field is based upon a com- 
bination of satellite data and surface gravimetry. Satellite solutions employing 
precise optical observations and Doppler data have been obtained by a number of 
investigators over the past several years. More recently, surface gravimetric 
data have been added to the solutions. Fields complete to degree and order 16 
have been obtained in this way. It is estimated that the satellite data contribute 
significantly to about degree and order 10. This corresponds to a half wave 
length of 18° , or about 2000 km. The uncertainty in the geoid height which is 
associated with the satellite contribution to this field is estimated to be on the 
order of 10 meters. This corresponds to an acceleration resolution of about 
2-1/2 milligals. 

It is planned to conduct a satellite-to-satellite tracking experiment between 
ATS-F and GEOS-C with a range accuracy of 2 meters and a range rate accuracy 
of 0.035 centimeters per second for a ten-second integration time. This experi- 
ment is now planned for 1974. It is anticipated that it will improve the spatial 
resolution of the satellite geoid by half an order of magnitude to about 6° . Longer 
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integration times should also permit a modest increase in the acceleration res- 
olution to about a couple of milligals. 

Satellite altimeter data will also be obtained by means of GEOS-C . An 
overall accuracy of five meters in altitude is the goal. The altimeter, per se, 
is expected to have an instrumental precision of about 2 meters, and an additional 
capability to observe with a precision of about 0. 2 meters for limited periods. 

The ability to interpret such measurements in terms of the sea surface topography 
and the gravity field is dependent upon the ability to represent the satellite's 
motion in altitude, to analyze the altimeter signal return in terms of sea state, 
and to model appropriately the departures of the sea level surface from the 
geoid due to tides, currents, the general circulation of the oceans, and atmos- 
pheric effects. Lasers with accuracies of about a decimeter are expected to be 
available in time to track GEOS-C. Three such lasers located at sites such as 
Antigua, Key West and the Canal Zone, for example, could yield satellite altitude 
determination uncertainties of the order of a meter over most of the tri ang ular 
region defined by these locations. It is estimated that satellite -to- satellite 
tracking can be used to extend the coverage over distances of oceanic scale with 
altitude uncertainties of less than 5 meters. Departures of the ocean surface 
from the geoid of up to one or two meters can be expected. Uncertainties 
associated with these increments are usually less than a meter over wide areas 
of the ocean. New information about the gravity field and the geoid with uncer- 
tainties in the range of a few meters should accordingly be derivable from the 
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GEOS-C altimeter data. This will represent an improvement in accuracy of a 
factor of two or so over many regions of the earth's surface. 

The detailed analysis of the satellite -to -satellite tracking and satellite altim- 
eter data will involve the simultaneous use of detailed mathematical models of 
the earth's gravitational field and figure, the tides, current systems such as the 
Gulf Stream, the general circulation of the oceans and, ultimately, factors 
associated with the ocean- atmosphere interface. Detailed analysis of the struc- 
ture and evolution of a tsunami will also involve ray tracing and modelling of the 
ocean's wave propagation characteristics and bottom topography. The develop- 
ment of mathematical models for geophysics promises to increase markedly in 
complexity as the satellite Earth Physics Program unfolds. 
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GRAVIMETRIC GEODESY AND SEA SURFACE TOPOGRAPHY STUDIES 
BY MEANS OF SATELLITE-TO-SATELLITE TRACKING 
AND SATELLITE ALTIMETRY 


I. INTRODUCTION 

The study of the earth's gravitational field by means of data from satellite- 
to-satellite tracking and satellite altimeter systems will take place in the context 
of the current knowledge of the subject. The consideration of these fundamentally 
new types of data will accordingly begin with a discussion of techniques which 
are already in use for probing the geopotential. The special nature of the ap- 
proaches offered by the new tracking capabilities will be the more clear as they 
are contrasted with the older, more conventional methods. 

The capability for studying sea surface topography opened up by the satellite 
altimeter is a two-edged one. The sea surface topography can be identified with 
the geoid, to within uncertainties associated with oceanographic effects and, con- 
versely, it can be viewed in terms of features of interest to the oceanographer, 
to within the errors associated with the gravity field representation. 

The actual positions of the sea sin-faces are affected by a number of ocean- 
ographic factors including, in addition to the geoid, the tides, the general circu- 
lation, currents and their meanders, sea-state, winds, storm surges and tsuna- 
mis. Hence, for example, in order to conduct the gravimetric investigations 
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properly, It is necessary to represent these oceanographic effects on the sea 
surface heights appropriately. Similarly, existing knowledge of the gravity field 
will play an important role in the investigation of oceanographic effects. Thus, 
gravimetric and oceanographic investigations relying upon altimeter data are 
interrelated. This point is discussed further later. 

The following discussion deals, in particular, with the potential contributions 
of the altimeter and satellite -to -satellite tracking systems associated with the 
GEOS-C spacecraft which is scheduled to be launched in mid-1974. The dis- 
cussion begins, actually, with consideration of studies which are either ongoing 
or getting underway and which should be completed by the time the data begin to 
flow from GEOS-C. These analyses will involve the use of the additional satellite 
tracking data, including especially the laser data, which will become available 
in the intervening years. This aspect of the investigation will involve the simul- 
taneous determination of gravimetric quantities and other environmental param- 
eters which contribute to the observational residuals. The aims of this part of 
the investigation will accordingly include the determination of tracking station 
locations, quantities such as the earth’s radius, and coefficients representing 
atmospheric drag and radiation pressure effects. Consideration is being given, 
too, to the possibility of equipping other spacecraft with the capability for 
satellite-to-satellite tracking through ATS-F. The potential contributions of 
such capabilities to the gravimetric geodesy investigations are also discussed. 
( 1 - 4 ). 
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II. THE EXPLOITATION OF EXISTING DATA IN PREPARATION FOR THE 


ANALYSIS OF ALTIMETER AND SATE LLITE - TO-SATE LLITE TRACK- 
ING DATA FROM GEOS-C 

A. General Considerations 

The basic strategy is to derive the greatest benefit from the altimeter and 
satellite-to-satellite tracking data from GEOS-C by analyzing it together with all 
of the related data which can contribute to the determination of the earth f s gravi- 
tational field. The determination of the characteristics of the gravitational field 
will be based upon the existing stores of geodetic data as well as the additional 
observational material which will become available in time for use in the analyses 
of the new types of data to be furnished by GEOS-C. The fundamental body of 
data upon which our current knowledge of the gravitational field is based consists 
of over a hundred thousand precisely reduced optical observations of some two 
dozen satellites, more than a hundred thousand electronic observations including 
both Doppler measures and range and range rate data, the beginnings of a supply 
of laser range measurements, as well as the gravimetric and survey results of 
classical geodesy. It is anticipated that, by the time the analysis of the new data 
types from GEOS-C commences, this store will be enriched through the addition 
of considerable numbers of laser observations obtained in programs following the 
ISAGEX campaign such as the San Andreas Fault Experiment (SAFE). 

It is contemplated that the interpretation of these collections of data will 
involve general determinations of the gravitational field, resonance determinations 
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of geopotential coefficients, and studies involving representations in terms 
of mass anomalies. 


Recent activity in the field of geopotential analysis at several institutions 
including Goddard is indicated in Tables I and II and discussed in references 28 
and 29 which include bibliographies. Results of comparisons of some of the 
geopotential models listed in Tables I and II are indicated in Table in and in 
Figure 1. These findings give some idea of the effective accuracies associated 
with some of the models. Improvements obtained by the modeling of resonant 
terms are also seen in Table in. 


Table I 

General Determinations of Geopotential Coefficients 


Author or 
Designation 

Year 

Reference 

Number 

Basis 

Field Charaet 

eristics 

Satellite Data 

Other 

Information 

Number of 
Satellites 

Tracking Systems 

Last term in 
Complete Portion 

Number of 
Coefficients 

1. NWL 5E-6 

1965 

1 

3 

Doppler 


7.6 

64 

2. APL 3,5 

1965 

2 

5 

Doppler 


8, 8 

84 

3. SAOM-1 

1966 

3 

16 

Cortical 


8, 8 

122 

4. Kaula K-8 

1966 

4 

12 

Optical, Doppler 


7. 5 

99 

5, Kaula C 

1966 

4 



Determinations 

7, 2 
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6, Rapp 

1967 

5 

16 

Optical 

Gravimetry 

14, 14 

219 

7, Kohnleln 

1967 

6 

16 

Optical 

Gravimetry 

15, 15 

250 

8. Kaula UCLA 

1967 

7 

9 

Optical, Doppler 

Gravimetry 

8, 8 


9. Rapp 

1968 

8 

16 

Optical 

Gravimetry 

14, 14 


10. SAO COSPAR 

1969 

9 

24 

Optical, Range, Range Rate 


14, 14 

280 

11. SAO B6.1 

1969 

10 

24 

Optical, Range, Range Rate 


16, 16 


12. SAO B 13.1 

1969 

11 

24 

Optical, Range, Range Rate 

Gravimetry 

16, 16 

314 

13. SAO SF 

1969 

12 

24 

Optical, Range, Range Rate 

Gravimetry 

16, 16 

316 

14, GSFC 1.70 C 

1970 

13 



Determinations 

15, 15 

249 






1-9, 11, 12 



15. SAO 69 (II) 

1970 

14 

21 

Optical, Range, Range Rate 

Gravimetry 

16, 16 

316 
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Table n 


Resonance Determinations of Geopotential Coefficients 



— I 

Reference 

Number 

Basis 

Field Characteristics 

Author or Designation 

Year 

Numbers and Types 
Satellites 

Reference 

Coefficients 

Coefficients 

Determined 

1. Wagner 

1967 

15 

3 24-hour 


<2, 2), (3,1), (3,3) 

2. Gaposchkin and Veis 

196? 

16 

3 12th order 

SAO M-l 

(13,12), (14,12), (15,12) 

3. Murphy and Victor 

1967 

17 

2 12-hour 


(2,2), (4,4) 

4* Yionoulis 

1968 

18 

3 13th order 


(13,13), (15,13), (17,13) 

5* Wagner 

1968 

19 

2 12-hour 


(3,2), (4,4) 

6. Murphy and Cole 

1968 

20 

1 12th order 


(14,12), (15,12) 

7. Wagner 

1968 

21 

3 24-hour 
2 12-hour 


(2.2) , (3,1), (3,2) 

(3.3) , (4,4) 

8. Douglas and Marsh 

1969 

22 

1 13th Order 

Yionoulis, 1968 

(14,13) 

9. Wagner 

1969 

23 

4 24-hour 
4 12-hour 

SAO B13.1 

(2,2), (3,2), (3,3), (4,4) 


Table D3 

Satellite Position Differences Associated with Various Gravity Models 


SAO M-l (modified*) 
vs. 

Position (meters) 

GEOS I 

GEOS n 

Radial 

Cross 

Track 

Along 

Track 

Total 

Radial 

Cross 

Track 

Along 

Track 

Total 

SAO M-l (unmodified) 

30.0 

17.1 

286.4 

288.5 

22.3 

10.3 

232.2 

234.2 

SAO COSPAR (no. llth) 

8.9 

12.8 

29.5 

33.4 

24.7 

18.2 

92.8 

97.8 

SAO 1969 

8.3 

13.5 

26.9 

31.2 

16.6 

19,3 

59.2 

64.4 

Kohnlein 

16.4 

16.1 

213.0 

214.2 

20.1 

16.9 

129.1 

131.7 

Kohnlein (modified) 

9.2 

10.9 

28.0 

31.4 

11.8 

15.1 

41.9 

46.1 

Rapp 

48.3 

29.5 

129.4 

141.2 

38,8 

37.7 

157.6 

166.6 

Rapp (modified) 

46.4 

33.2 

99.9 

115.0 

36,2 

39.1 

84.8 

100,1 

APL 3.5 

46.1 

46.8 

175.7 

187.1 

71.8 

55,2 

674.3 

680.4 

APL 3.5 (modified) 

42,5 

41.6 

90.1 

107.9 

34,7 

45.4 

88.1 

105.2 

NWL 5E-6 

16.3 

16.6 

204.0 

205.3 

46.4 

80.9 

374.6 

386.1 

NWL 5E-6 (modified) 

16.7 

12.9 

49.1 

53.4 

26.3 

80.0 

82.9 

118,2 

Kaula 

32.5 

42.2 

114.1 

125.9 

47.6 

43.5 

232.8 

241.4 

Kaula (modified) 

32.1 

42.5 

110.2 

122.4 

48.7 

42.0 

140.5 

154.6 


*Gaposchkin & Veis (1967) 12th order terms for GEOS-1 and Yionoulis (1968) 
and Douglas & Marsh (1968) 13th order terms For GEOS-1 1. 
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The program for the study of earth's gravitational field at the Goddard Space 
Flight Center has several major aspects. These come together in a program for 
a general determination of both the zonal and tesseral harmonics of the geopo- 
tential, the locations of tracking stations, and fundamental parameters such as 
the earth's radius. These quantities are determined simultaneously in this portion 
of the program. In addition, the overall Goddard program includes specialized 
research efforts directed toward the determination of zonal harmonic coefficients, 
resonant tesseral terms both individually and in sets, tracking stations locations 
and the geoid. 

B. The Goddard Earth Models 

The program for the general determination of the gravitational field has 
thus far generated four earth models, (30) These are referred to as elements 
of the Goddard Earth Model (GEM) series, individual solutions being designated 
GEM-1, GEM-2, etc. A broad view of the four GEM results obtained so far is 
seen in Table IV. More details about the GEM-1 and GEM-2 solutions are sum- 
marized in Table V. The body of optical satellite tracking data which constitute 
the principal basis of these results is described in the listings of Table VI. The 
surface gravity data which permit the filling in of fine structure detail in many 
areas are portrayed summarily in Figure 2. Information related to the un- 
certainties associated with these two types of data is seen in Figure 3. The 
relative weighting of the satellite and gravimetric data in the combined solution 
was based largely upon this type of information. 
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Table IV 


Description of Goddard Earth Models (GEM) Geopotential Solutions 


Solution 

Spher. Harm. 
Field 

Satellite Data 

Gravimetry Data 

GEM 1 

12 x 12 
Satellite 

120,000 Optical Obs. 
on 25 Satellites 


GEM 2 

16 x 16 
Combined 

GEM 1 Data 

1707 5° x 5° Mean 
Gravity Anomalies 

GEM 3 

12 x 12 
Satellite 

300,000 Optical and 
Electronic Obs. 
on 27 Satellites 


GEM 4 

16 x 16 
Combined 

GEM 3 Data 

1707 5° x 5° Mean 
Gravity Anomalies 


• 61 center of mass tracking station locations 

• Two low inclination satellites SAS 3°, PEOLE 15° 25 satellites in previous 
solution all have inclination >28° 

• 5° x 5° mean gravity anomalies based upon 
19,000 1° x 1° Anomalies from ACIC 

2,000 1° x 1° Anomalies from other sources 
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Table V 


Summary of Solution, Data, and Results 


A. SOLUTION 

Geopotential solution in spherical harmonic series (degree n, order m). 


Type Solution 

Complete 
(n x m) 

Zonals 

Complete 

Resonant Coefficients 

Satellite (GEM 1) 

12 x 12 

21 

m = 12, 13, 14, n = m to 22 
m = 9, n = 9 to 16 

Satellite & 



m = 11, n = 11, 12, 14 

Gravimetry 
(GEM 2) 

16 x 16 

21 

Complete as above 


46 Center of mass tracking station locations (see tracking systems 
below) 


B. DATA 

Satellite Data 

Tracking Systems 25 Satellites (Data Processed) 

300 weekly arcs of optical data 
15 weekly arcs of electronic/laser data 
66 one/two-day arcs of GEOS-I and n 
flashing light data 
20 weekly arcs of Minitrack data 

Gravimetry Data 

19,000 1° x 1° mean gravity anomalies of ACIC data, and 
2,000 1° x 1° mean gravity anomalies from other sources, used to 
form 

1,705 4° x 5° means, providing 70% world coverage 
C- RESULTS 

1. Station heights indicate a e = 6378145 meters. 

2. Mean value of gravity (Ag, = 3.28 mgal) indicates a e = 6378142 (m). 

3. Solution for 327 zonal and tesseral coefficients including Ag 0 
from satellite and gravimetry data. 

4. Solution using optical and electronic tracking data for location 
of 46 tracking stations . 


13 Baker -Nunn Cameras 
23 MOTS Cameras 
3 GRARR (range/range -rate) 
2 LASERS 
5 DOPPLER 
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Table VI 


Satellite Data 


300 WaaLly Saving Arc* Qpl. Data (Primarily SAP Bohw-Nunn> 


SotallTtc 

Mama 

A 

(Mafari) 

E 

1 

<0*g} 

ParlQ&a 

Haight 

IM 

Parted 

(Rav./Doy) 

No, 

Arct 

No. 

Obi. 

Avg. No. 
Oh* ./Are 

Avg. RMS 
(Wflightotfl 

TELSTAR-I 

9469530.1 

0.2421 

44.79 

951.3 

9.13 

14 

1944 

121 

2.70 

TIROS- 9 

BG20761 .2 

0.1147 

96.42 

704.7 

12.09 

14 

1525* 

109 

1.14 

GEOS-l 

B067353.6 

0.0725 

59.37 

I107.J 

11.96 

35 

45555** 

1301 

1 .26 

5ECOR-5 

8154849.9 

o.oeoi 

69.23 

1140. 1 

11,79 

4 

290 

72 

2.38 

OV1-2 

B314700 .2 

0.1635 

144,27 

414.6 

11.45 

4 

91 D 


1.93 

A LOU-2 


■’Rkyf 


dZ&9 

11.91 

6 

590* 

. 98 

0.B9 

ECHO-IRB 


WMM 


WPiM 

12.20 

18 

2240 

124 

2.24 

Dl-D 




599.0 

13.07 

9 

4386 

709 

1.64 

BE-C 




941 .9 

13.34 

22 

4947 

224 

1 .59 

Dl-C 



■£2121 

9E229 

13.79 

4 

902 

225 

2.53 

ANN A- IB 





i3, as 

40 

4183 

104 

1.5) 

GEOS-2 

KSEZZQi 

bR'> r ':M 

105.79 

1114.2 

12.62 

24 

25315** 

1034 

1.75 

OSCAR-7 

(ZZZSEI 


(□23 

847.7 

13.63 

^ 1 

1780 

445 

2.34 

S6N-2 

(222S9 


■o 

1062.5 

13.47 

5 1 

1SS 

71 

5.17 

COURIER-i B, 

(22H22M 


K23 

968.5 

13.44 

12 : 

3375 

261 

1.66 

GRS 


RO 

mmm 

K(9 

14.13 

5 

369 

73 

3.04 

TRANSITS 



(EEM 

■£23 

13.86 

14 

131* 

94 

5.92 

6E-B 



WESSM 

9E2E9 

13.74 

4 

469 

117 

1.87 

OGO-2 

IEIESEEI 


87.37 

424.8 

13.79 

7 

461 

65 

3.47 

INJUN-1 



(29 

WESEM 

13.68 

9 

768 

! 65 

2.15 

agena-r 


■.'JT'I’IL'M 



13.93 

7 

1003 

143 

2.84 

Ml 0 AS -4 


■ WlUB 

(C32B 

^n<:w 

0.69 

20 

14879 

1 743 

1.20 

VANGUARD-2 


0.1645 

32.89 

5*4.7 

11.47 

IL 

379 

34 

1.13 

VANGUARD-2S 

6309120.5 

0.1448 

32.67 

542.2 

11.46 

5 

615 

123 

2.29 

VANGUARD-3S 

SSI 1504 .6 

0.1904 

33.35 

517.9 

11.06 

15 

996 

64 

2.09 

Summery Tote It | 

314 

121554 

3S6 

1,4* 


“MOTS 

»aocnc*n 


GIOSII. II 


The Goddard Earth Models were determined by means of the method of 
numerical integration. The program system used is based on the one described 
by Velez et. al. (31) The special perturbation method employed here is funda- 
mentally different from the general perturbation technique used by SAO in gen- 
erating the 1969 Smithsonian Standard Earth (II)* 18 \ Sharper insights are pro- 
vided by the latter, analytical, approach. Increased accuracy is usually more 
readily achievable by means of the former, numerical method. 

The actual values determined in the GEM-4 solution are seen in Table VII 
which contains a listing of both zonal and tesseral geopotential coefficients, and 
Table Vin which lists the tracking station locations. The geoid heights corres- 
ponding to the values of Table vn are portrayed in Figure 4. The corresponding 
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USED FOR RELATIVE WEIGHTING 
IN COMBINED SOLUTION 



DEGREE n 

Figure 3. Average Standard Deviation of Coefficients per Degree n of Satel- 
lite and Gravimetry Solutions 
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Table VII 

GSFC Geopotential Solutions (Normalized Coefficients x 10 8 ) 1 



L 

M 

GEM 4 

GEM 3 


L 

M 

GEM 4 

GEM 3 


L l 

N 

GEM4 

GEM 3 


L 

m 

GEM 4 

GEM 3 


L 

M 

GEM 4 

G*M) 

c 

2 

r 

>484. 189.: 

-484. 1718 

c 

21 

0 

-C .Co 76 

-C.0CT5 

C 

4 

2 

0*3511 

01.3942 

C 

9 

3 

-0*1760 

-5. 1403 

c 

15 

4 

0*0099 

0.0 

s 

2 

C». 

<)*>) 

0.0 

s 

21 

0 
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GEM 4 Station Coordinate Solution 
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contours derived from the SAO 1969 (II) solution are also shown there. Some 
feel for the accuracy of the GEM-4 solution can be had by inspecting Figure 5. 

The closures of the tracking station heights have a scatter of the order of ten 
meters about a mean which corresponds to the value of 6378145 meters for the 
earth’s equatorial radius. 

C. Zonal Harmonics 

Zonal harmonic coefficients obtained simultaneously with the tesseral har- 
monic coefficients in the GEM solutions are based on arcs a week in length. In 
the past zonal harmonics have been derived separately and, in fact, are often 
used as fixed reference quantities in the determination of the tesseral harmonic 
coefficients (32 » l3) . Zonal Harmonic coefficients have secular and long period 
effects on all the elements except the semi -major axis. Collections of elements 
spanning at least a complete revolution of the perigee accordingly provide a good 
basis for zonal harmonic solutions. The earth’s pear- shape was discovered in 
this way through examination of the eccentricity perturbations of the Vanguard I 
satellite. Kozai and King-Hele have extended this process using theoretical 
approaches which are primarily analytical in nature. Cazenave, et. al. , using 
a similar approach, derived an extension of Kozai 's field by adding data for three 
low inclination satellites. ( 3 3) Wagner has recently derived a solution using a 
method of numerical integration of the orbital elements. (34) It is of special 
interest to compare Wagner’s results derived by the classical approach, which 
relies on long arcs whose lengths spans a full revolution of perigee, with the 
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GEM-2 results which are based on arcs an order of magnitude shorter, i.e. , 
about a week in length. This is done in Figures 7 and 8. The field derived by 
Cazenave, et. al. , (33) is also indicated in the comparisons in these figures 
where it is designated FR 71. The basis for each of these three solutions is 
tabulated in Table IX and depicted graphically in Figure 6. Results and estimates 
of uncertainties are indicated in Table X and XI. 

The SAO field listed in Table X was obtained using orbits and a method 
similar to those employed in deriving the FR 71 set of coefficients. Compari- 
sons of Wagner's field with GEM-2 and FR 71 results in terms of secular rates 
and eccentricity oscillation amplitudes are seen in Figures 8 and 8. Finally, 
zonal geoid profile information is given in Figures 9 and 10. 

D. Resonant Tes serai Harmonics 

Resonant tesseral harmonics also can be studied very effectively through the 
analysis of longer arcs. 

Some of the results obtained by Wagner, Murphy, and Victor are contained 
in references 19, 21, 23-25, 27, and 34A, summarized in references 28 and 29, 
and indicated in Figures 11 through 14 and Tables XII through XIV. 

The possibilities for obtaining an improved geopotential solution through the 
more extensive use of synchronous and semi- synchronous satellite data are in- 
dicated in Table XIV where it is seen that relatively small changes in the coef- 
ficients result in a significant improvement in the fits to the orbital data. A 
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Satellites Used in Recent Zonal Solutions 
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Figure 6. Satellites Used in Present Solution 
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Table X 


Zonal Coefficients from Recent Solutions 
(unnormalized in units of 10 -6 except as noted) 
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Table XI 

Total Weighted RMS Residuals in Road Solutions 
for Long Satellite Arcs 
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SAS 1 

EXPLORER I 
ANNA IB 


TELSTAR I 

GEOS 2 

EGRS 3 


PEGASUS 3 


FR-1 


OSO 3 


EXPLORER 27 


ESSA 1 


TIROS 5 


TIROS 9 


ISIS 1 


AVG. RMS IN 21 ARCS 


WITH FIELD USED 


PRESENT 
SOLUTION 
{WAG 72} 


1.52 


CAO 2 
(a = 1.12 

e - 0.001) 

1965-34A 


(a = 1.51 

e - 0.05) 

ALOUETTE 2 

(a = 1.27 

e = 0.15) 

COSMOS 382 

(a = 1.60 

e = 0.14) 

DIAL 


(a = 1.15 

e = 0.09) 

EGRS 5 {SECOR 5} 

(a = 1.28 

e = 0.08) 



1.65 


6 SATELLITES NOT USED IN PRESENT SOLUTION 
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Figure 9. Zonal Geoid Profile for WAG 72 
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Figure 10. Differences in Zonal Geoid Profiles 
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Figure 11. Semi-Major Axis of 1965 30A [Molniya] vs. Time 


substantial improvement is obtained in residuals for synchronous and semi- 


synchronous satellites by solving for five pairs of tesseral harmonic coefficients 


in the context of a recent SAO general field. Improved values for low order 


harmonics which represent these features, and for the general field as a whole, 


can accordingly be expected to result from the inclusion of these added data in 


subsequent solutions in the GEM series. 


Higher order resonances have been studied by somewhat different methods* 
Use has been made, for example, of general perturbation approaches such as 
those discussed in references 21, 24, 28 and 25 through 37. 

An example of results of such a determination of resonant harmonic coeffic- 
ients based on Minitrack data is indicated in Figure 15. The influence of this 
satellite is also reflected in overall general solutions in the GEM series. 
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LONCmiDf DRIFT ACCELERATION: 10 * * RADIANS /SI D£ft£AL DAY* 
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Figure 12. Drift Acceleration in Syncom 3 Arcs 6, 7, 10, 11 and 13 
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Figure 13. Drift Acceleration and Estimate of Acceleration 
Uncertainties on 24 Hour Satellites 
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Figure 14. Observation Residuals for 12- and 24-Hour Satellite Accelerations 

from Recent Geoids 
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Table XH 


Resonant Longitude Gravity Coefficients* from Recent Studies 


Geoid 

C 22 
< io - 6 ) 

s 23 

( 10 *) 

Oil 
( 10 '* ) 

s,. 

( 10 -*) 

c 32 
( 10 -®) 

s 32 

( 10 -*) 

( 10 -*) 

s 33 

( 10 -*) 

( 10 -*) 

S 4J 

( 10 -*) 

C 44 

( 10 -*> 

S 44 

( 10 -*) 

Kohnlein 

2.38 

- 1.35 

1.71 

0.23 

0.84 

- 0.51 

0.66 

1.43 

0.35 

0.48 

0.04 

0.30 

Kaula [1966] 

2.42 

- 1.36 

1.79 

0.18 

0.78 

0.75 

0.57 

1.42 

0.30 

0.60 

- 0.06 

0,32 

Wagner 

2.40 ± .03 

- 1.43 ± .03 

- 0.42 

-1.58 

0.69 ± .20 

- 0.53 ± .20 

0.16 

1.10 

0.0 

0.0 

0.02 

0.70 


*Ful!y normalized coefficients. 


Resonances have also been determined by means of a numerical integration 
approach. Results of this kind are indicated in Figures 16 and 17 in reference 
26. 

E. Tracking Station Locations 

Dynamical determinations of tracking station locations were carried out in 
the Goddard Space Flight Center prior to the derivation of the GEM 
solutions/ 38 * 41 1 Some indications of the coverage provided by the thirty thousand 
observations employed in a recent determination is indicated in Table XV and 
XVI and Figure 18. Results are presented in reference 41. 

F. Geoid Studies 

A detailed gravimetric geoid from North America to Eurasia has been de- 
rived by Vincent, Strange and Marsh by combining information resulting from 
both satellites and surface gravity studies. 1421 A comparison of this result with 
the geoid based on the SAO 69 (II) field is seen in Figure 19. 
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Table XIII 

24 Hour Satellite Arcs in 1970 Resonant Geopotential Solutions 


SATELLITE 

ARC 

NUMBER 
OF KEPLER 
ELEMENT 
SETS 
USED 

SPAN OF 
DATA 
(modified 
Julian 
Days) 

GEOGRAPHIC 
LONGITUDE 
SPAN IN 
ARC 

(degrees) 

ORBIT 

INCLINATION 

(degrees) 

ESTIMATE OF 
QUALITY OF 

RMS RESIDUALS IN 
LONGITUDE FROM 
COMBINED ARC SOLUTIONS (degrees) 

t) E 1 ERMlNAT ION 
OF THE GEOGRAPHIC 
LONGITUDE. </ (X), 

IN ARC, FROM 
INDIVIDUAL ARC 
ANALYSES 
(degrees) 

WITH 24 

HOUR SATELLITES 
ALONE, 

SOLVING FREELY 
FOR 6 RESONANT 
COEFFICIENTS 

WITH 24 HOUR 
SATELLITES 
COMBINED WITH 
SAD 1969 COSPAR 
DATA, !N CONSTRAINED 
SOLUTION FOR 10 
RESONANT COEFFICIENTS 

WITH SAO 
1969 COSPAR 
FIELD 1 

ALONE 

SYNCOM 2 

1 

16 

3826348351' 

302-305 

33.0 

0.025 

0.013 

0.013 

0.013 

SYNCOM 2 

2 

11 

333G1-3B443 

236 301 

32.8 

0.025 

0.034 

0.034 

0.035 

SYNCOM 2 

4 

9 

38510-38570 

19G-243 

32.5 

0.020 

0.015 

0.014 

0.014 

SYNCOM 2 

5 

24 

38580-30807 

72-189 

32.3 

0.035 

0.030 

0.033 

0.035 

SYNCOM 2 

8 

15 

3881G-3391& 

65-68 

31.8 

0.020 

0.018 

0.018 

0.018 

SYNCOM 2 

DoD 

42 

38816-40104 

65-86 

29-32.0 

0.040 

0.037 

0.038 

0.103 

SYNCOM 3 

6 

9 

38690-38750 

178-180 

0.0 

0.015 

0.014 

0.015 

0.013 

SYNCOM 3 

7 

10 

38775 38335 

174 181 

0.0 

0.055 

0.054 

0.055 

0.052 

SYNCOM 3 

11 

16 

30070-39262 

165-172 

0.5 

0.020 

0.012 

0.014 

0.013 

SYNCOM 3 

13 

4 

39385 39491 

1G0-1G1 

1.3 

0.015 

G.G05 

0.006 

0.010 

SYNCOM 3 

14 

15 

3DGG3-4Q175 

146-160 

23.0 

0.045 

0.001 

0.061 

0.110 

ATS 1 

1 

12 

40204 40248 

210 211 

1.2 

0.005 

0.001 

0.001 

0.001 

ATS 3 

1 

12 

40197-40241 

288 312 

0.6 

0.005 

0.002 

0.003 

0.003 

ATS 3 

2 

13 

40267-40337 

237 288 

0.3 

0.005 

n.oo4 

0.003 

0.002 

ATS 3 

3 

10 

40524-40577 

313-315 

0.3 

0.005 

0.002 

0.003 

0.003 

INTELSAT 2*F3 

1 

12 

39G07-39905 

349-352 

1.0 

0.025 

0.028 

0.028 

0.030 

INTELSAT 2-F3 

2 

10 

4040G40G42 

34G-349 

1.0 

0.020 

0.011 

0.012 

0.015 

INTELSAT 2-F4 

1 

8 

40323-40GQ8 

179 194 

1.0 

0.040 

0.043 

0.046 

0.062 

EARLY EURO 

1 

11 

38897*39080 

324 332 

0.5 

0.025 

0.025 

0:026 

0.026 

EARLY BIRD 

2 

6 

3909G-39Z19 

321-331 

0.8 

0.010 

0.007 

0.009 

0.009 

ATS 5 

1 

14 

40476 40556 

253-256 

2.6 

0.005 

0.004 

O.OOG 

0.009 


NASA GSFC-T 8. OS 




Table XIV 


Resonant Geopotential Coefficients for 24 Hour Satellites and RMS Fits to 21 Arcs 
: UNNORMALIZED COEFFICIENTS IN UNITS OF 10‘ 6 : 


Field 

Overall 
Weighted 
RMS 
Residual 
in 21 
Arc 

Solution 

2,2 

3,1 

3,3 

i 

4,2 

4,4 

C 

S 

C 

S 

C 

S 

C 

S 

C 

S 

SAO STANDARD 
EARTH 11 (1970) 

2 Ad 

1.558 

-0.881 

2.12S 

0.281 

0.096 

0.199 

0.074 

0.158 

-0.0017 

0.0072 

SAO CO SPAR 
(1909) 

1.38 

1.560 

-0.890 

2.040 

0.2G2 

0.096 

0.198 

0.073 

0.148 

-0.0028 

0.0078 

WAGNER (UN- 

j CONSTRAINED) 

! 2,2— 3,1— 3,3 (1970) 

1 

0.849 

1.56S 

-0.907 

1 

! 

1.687 

0.483 

0.103 

0.204 

0.074 

0.158 

-0.0017 

0.0072 

WAGNER— SAO 
COSPAR (CON- 
STRAINED) (1970) 

0.S9G 

1.570 

-0.908 

2.029 

0.267 

0.098 

0.205 

0.075 

0.150 

-0.0029 

0.0079 



SEMIMAJOR AXIS (ka) 
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808.623 


808620 
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Figure 15* Semimajor Axis of Tiros EX 


X 10H 

£ ZERO LINE 


♦ ROSRAN 

* ALASKA 

O RANGE xCARVON 

• R RATE + MADGAR 

AMPLITUDE = 23 MILLISECONDS 


MAU. ftW.T IM 

MOSIOM * TWjKTfC*V ANAiVSft DNftKX 
MUNCH |»Q 0AT1 i/Tfl 
** AUftftl NO. 1*4 


Figure 16* Apparent GEOS-II GRARR Timing Errors SAO Ml Gravity 
+ APL 13th-Order Terms 4/29/68-5/4/68 
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Figure 17. Apparent Rosman Range Timing Errors with SAO Ml Gravity, 
APL 13th-Order Coefficients and C 14i 13 = 0.57 x 10 -21 , 

S 14 13 = 6.5 x 1CT 21 , 4/29/68 to 5/4/68 


HI. THE ANALYSIS OF GEOS-C DATA 

A. Satellite -to -Satellite Tracking Data 

1. The Satellite -to -Satellite Tracidng System 

The satellite-to- satellite tracking experiment involving a geodetic 
satellite calls for the tracking of GEOS-C via ATS-F and, later ATS-G. The 
quantities to be observed are the range from a ground station through ATS-F to 
GEOS-C, and the corresponding range difference over a small time interval 
which is often referred to as range rate. The vise of these data involves the 
analysis of the motions of both the GEOS-C and the ATS-F spacecraft. It will be 
possible to make independent measurements of the range and range rate from the 
ground tracking station to ATS-F using the basic ATS tracking system. It is 
anticipated that the ATS-G spacecraft may also be fitted with laser corner re- 
flectors, and hence that this accurate tracking system can also be used to obtain 
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Table XV 


.Number of Optical Observations Per Station Used in Dynamical Estimation 


Station 

Observations 

Blossom Point, Md. 

(1021) 

650 

Ft. Myers, Fla. 

(1022) 

680 

Woomera, Australia 

(1024) 

152 

Santiago, Chile 

(1028) 

264 

Mojave, Cal. 

(1030) 

932 

Johannesburg, Un. of So. Afr, 

(1031) 

623 

St. Johns, Newfoundland 

(1032) 

179 

E. Grand Forks, Minn. 

(1034) 

542 

Winkfleld, England 

(1035) 

400 

Rosman, N. C. 

(1037, 1042) 

1015 

Orroral, Australia 

(1038) 

482 

Tananarive, Madagascar 

(1043) 

339 

Univ. of N. Dakota, N. D. 

(7034) 

548 

Edinburg, Texas 

(7036) 

1109 

Columbia, Mo. 

(7037) 

1540 

Bermuda 

(7039) 

448 

San Juan, P. R. 

(7040) 

475 

Greenbelt, Md. 

(7043) 

158 

Denver, Colo. 

(7045) 

573 

Jupiter, Fla. 

(7072) 

516 

Sudbury, Ontario 

(7075) 

699 

Kingston, Jamaica 

(7076) 

388 
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Table XVI 

GRARR and Laser Arcs Used in Solutions 






H?6fl 

ARCS 





i 

i 


I'LASKH 


MAUGAX 


LOSHAN 

WALLAS 

GODLA0 I 

Dale 

' No. of Os. 

Nn. of Ot a . 


Vo. of ! 

b*/. 

No. of 

PSB.ieB 

N of 
Ob*. 

F< an^e 

No, of 


■I —c" 

Rargrc V** so* 

Kang© 

i U* p ”**' 

15 ■nr 7 f j 

i-VC 

l»l!J8CS 

i«*. Pa,t€ ' f 


215 H9 3 




180 

IHO 

3 

495 

2 


4/26-27/68 

136 84 3 




313 

313 

4 

U9 

1 


5/7-8/G8 

172 172 2 




102 

102 

1 

460 

2 


5/21-22/68 

132 132 2 




157 

ICO 

4 

402 

2 


e/9-1 Q/68 

262 252 3 




2T7 

27 7 

4 




6/11-1 2/08 

162 182 3 




456 

456 

S 

508 

3 


6/14-15/68 

252 252 4 




1M 

ID? 

3 




6/1G-17/C8 

136 136 2 




193 

194 

3 




6/21-22/06 

271 271 5 




2C3 

279 

4 

422 

3 


6/23 -*4/68 
9/24-20/68 

285 285 5 

111 

112 

2 

229 

229 

4 



69 1 

0/27-28/68 

193 133 3 


51 

2 






38 1 

10/4-5/68 

200 200 3 

105 

1 07 

2 






73 I 

10/6-7/08 

271 271 4 

1 86 

189 

4 







10/8-9/08 

67 67 1 

117 

169 

4 






3!8 2 

10/21-22/68 

132 132 2 









212 2 

1 C/23- 24/60 

202 202 3 









437 2 

TOTALS 

3108 2990 48 

579 

653 

14 

2375 2393 

34 

2530 

13 

1247 0 





1069 

ARCS 





i 


U LAStCR » 



CARTON 



CKI.fl/l 8 

GOUT AS' | 

Date 

No- ol Obfl . 


No. of Obfl. 

No. of 
PfcftfiCft 


No. of 
Ob*. 

Kin?e 

No. of 
PiUiflc* 

No. of 
Obs. 

fl ange 

No, of 
Paflfion 

1*0. of 

flange Piser* 

flange 

R /trifle 

n.-tr 


3/2-3/69 

127 127 

5 

150 

138 

4 


9S 

2 

64 

1 

3/5-6/69 

124 124 

3 

92 

123 

2 


227 

J 

3C6 

4 

3/11-12/C9 

99 99 

3 

232 

23? 

4 


190 

3 

106 

2 

3/13-14/09 

150 150 

3 

314 

379 

6 


190 

3 

32 

1 

3/17-1 8/G9 

1% 196 

b 

232 

237 

S 


214 

4 

234 

3 

3/29-30/69 

92 92 

5 

170 

170 

S 






3/31-4/1/69 

164 . Ifi3 

4 

192 

200 

1 


101 

2 


t 

4/6-9/69 

75 75 

2 

167 

176 

4 


14C 

2 

91 

3 

4/10-1 1/G9 

172 172 

5 

90 

113 

2 


321 

4 

73 

1 

4/14-15/G9 

123 123 

3 

99 

125 

3 


271 

4 



4/24-25/69 

139 159 

5 

ISO 

179 

3 


251 

£ 

155 

1 

5/5-6/69 

163 163 

4 

219 

280 

4 


216 

4 



TOTALS 

1684 1683 

47 

2115 

2414 

43 


2235 

33 

1?27 

16 

•Station coordinate* held fixed. 











SUMMARY 


GRARR 

No. of Qbs. 

CRARR 

No. ofjObs. 

range 

9302 

Laser 


range rule 

10133 

range 

7205 

number of pasnes 

185 

number of passes 

71 


NAME 

LOCATION 

INSTRUMENT 

ULA5KP 

Kalrbanks, Alaska 

GTl/.RR 

MADGAH 

TanariATivt , 

CRARJl 

ROSUAN 

Ac*in?.r, N.C. 

GRARH 

CATiVON 

Carnarvon, Austral'* 

CRARJR 

WALLAS 

Wallop* 1*1 pud. Virginia 

Loser 

GOD LAS 

GrrtniiHi, Mart land 

l.JLRer 

CKMLAS 

Carnarvon, Australia 

Larn? 
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Figure 18* Subsatellite Plot of GEOS-I and GEOS-II Passes, Edinburg, Texas 
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Figure 19. Detailed Gravimetric Geoid and SAO ’69 Geoid 
(Gaposchkin and Lambeck, 1970) 
at 38° N Latitude 

The tracking of GEOS-C from ATS will be analogous in many ways to the 
tracking of lunar satellites from stations on the earth. The data from lunar 
orbiters were originally analyzed in conventional ways in terms of spherical 
harmonic coefficients. It was realized after a time that many of the lunar satel- 
lite range rate tracking data residual anomalies correspond to mascons or sur- 
face density distributions. This interpretational approach is perhaps one of the 
most striking features of recent research relating to gravitational fields. (Cf. 
Figure 20). Both the traditional and the new approaches have been used at 
Goddard as is pointed out in reference 43 and 44. 
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Figure 20. Doppler residuals from Lunar orbits of Apollo VIII on Decem- 
ber 24, 1968. The upper and lower portions of the Fig. depict residuals 
obtained in the neighborhood of Sinus Aestuum during the first and seventh 
revolutions, respectively. The residuals are plotted vs. position along the 
retrograde orbit. The positions in the two revolutions are indicated by 
surface tracks which are correspondingly marked in the central portion of 
the Fig. The Apollo Lunar landing sites are also shown. 
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The tracking of a low altitude satellite such as GEOS-C from ATS-F is ex- 
pected to yield information about the gravitational field which is more detailed in 
terms of space resolution than that obtainable from the existing sets of satellite 
data. The degree to which mean gravity anomalies in adjacent squares on the 
earth's surface can be separated on the basis of satellite-to-satellite tracking 
data of the ATS/GEOS-C type is a function of the altitude of the low satellite. 

The dimension of such a square are roughly comparable to the altitude of the 
low satellite. This problem has been studied in some detail by Schwartz* 45 ^ . 
Figure 21 depicts the relationship he found between the satellite altitude and the 
size of the square in which gravity anomalies can be identified by means of 
satellite-to-satellite tracking techniques. Extrapolation of this result to the 
general altitude region now planned for GEOS-C, i.e. , the neighborhood of 800 
to 1000 kilometers, leads to an estimate of the order of six degrees for the 
spatial resolution capability of the GEOS-C /ATS-F experiment. 

i 

The magnitudes of the effects are indicated in Figure 22 for a satellite in a 
700 km altitude orbit. It is presently estimated that the ATS-F/GEOS-C tracking 
system will have an accuracy of about 0. 35 mm/sec for a twenty second integra- 
tion time. The corresponding resolution capability in terms of a mean anomaly 
in a five-degree square for a satellite in a 700 kilometer altitude orbit is about 
4.4 milligals. The implications of this for GEOS-C will be discussed further in 
Section IV. 
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2. GEOS-C Orbit Selection Considerations 

The following discussion is intended to indicate the different kinds 
of scientific objectives which can be met with a particular selection of a GEOS-C 
orbit. 

It is assumed for the purpose of the case considered here that the GEOS-C 
orbit will have an inclination of about 115° , i.e. , 65° retrograde and a mean 
altitude in the neighborhood of 800 to 1000 kilometers. The eccentricity will be 
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Figure 22. Sensitivities of p in thousandths of a millimeter per second to in milligals. One satellite in a 
700 km high orbit is tracked by a very high geostationary satellite. E denotes the position of the low satellite 
at the initial epoch, and P denotes its position at the epoch for which the partial derivatives are evaluated. 



taken to be 0.006 or less, a range which will simplify the altimeter design. One 
is at liberty to adjust the mean altitude or the period within certain limits in the 
attempt to achieve as many of the scientific goals as possible. The values finally 
selected will be governed partly by considerations brought out in the following 
discussions. They may be somewhat different than those indicated here, however 
it is not anticipated that the principles underlying the present discussion will be 
fundamentally changed. 

The spatial resolution obtainable with the satellite-to- satellite tracking data 
in an experiment of the GEOS-C/ATS-F type is about 5. 7 degrees for the case 
considered. Starting from this point, one can work out the following orbit param- 
eters which will accommodate the objectives of a gravimetric geodesy investi- 
gation from the standpoints of the gathering of both the satellite-to-satellite 
tracking data dn altimeter tracking data. It will be seen later that a number of 
important oceanographic objectives can also be met with the orbit selection which 
is discussed here. 

An orbit having a 65° retrograde inclination, a mean altitude of approximately 
1000 km and a corresponding period of about 105 minutes will move westward in 
longitude by about 26.2° each revolution. This orbit has ground tracking equator 
crossings which shift to the west about 6.3 degrees per day. The resulting 
separation between the orbit traces, which is smaller by the factor sin i, is 
5.7° which corresponds to the value indicated above. The desired 5.7° survey 
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pattern is thus approximately traversed in fifty-five revolutions of GEOS-C, on 
a bit more than four days. It includes a complete set of tracks, crossing the 
equator in both the ascending and descending directions separated by about 5.7 
degrees. 

At the beginning of the fifth day, the longitude of the ascending node of the 
56th revolution is found to be displaced from that of the first revolution by about 
a degree, i.e. , about a sixth of the daily shift. (Cf. Figure 23). This can be 
thought of as the beginning of the phase in which the fine, or one-degree, pattern 
is laid down. The process continues for a total of 344 revolutions, at the end of 
which the ground track repeats, i.e. , the ground tracks of the first and the 
345th revolutions would coincide. This, in fact, provides a precise defining 
criterion for the orbit in the context of the chosen inclination, i.e. , 65° retro- 
grade, and eccentricity, namely, about 0.006. Thus, in something over 25 days, 
such a GEOS-C orbit would provide for an altimeter survey with orbits separated 
by slightly less than a degree. Since the altimeter cannot operate continuously, 
due to power limitations, an actual survey of this type would take much longer, 
perhaps on the order of a year or more. (46) 

Clearly other strategies are possible, e.g., by selecting patterns which 
would give spacings of about 3°, 1.5°, etc. The example sketched here will 
suffice for the purposes of the present discussions, however. Resonances may 
be associated with some of these choices. A preliminary look at this point 
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Figure 23. Equator Crossing Pattern for a Typical GEOS-C Orbit 
indicates that these will not be unduly severe, however. The degree to which such 
a specific fine resolution can be realized depends upon the launch vehicle 
capabilities. 
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3. Satellite-to-Satelllte Tracking Data Requirements 
a. Gravimetric Geodesy Analyses 

This information about the orbit furnishes the basis for esti- 
mating typical data taking requirements. The 6.3 degree separation at the 
equator implies that some 57 passes will be trackable from ATS while it is in 
a fixed position. 

Half of these correspond to ascending nodes, half to descending nodes. Each 
pass will be about 58 minutes long. Hence, a total of some 56 hours of satellite- 
to-satellite tracking would complete such a six-degree survey. Redundant in- 
formation would be obtained near the maximum latitudes since the tracks are 
less than 6 degrees apart. In an initial experiment, however, it would probably 
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be wise, even so, to plan on continuous satellite-to-satellite tracking which 
should aid considerably in the interpretation of the data. 

ATS-F will be at 94° W longitude at first, and then at 35° E longitude. Thus 
some 112 hours in all would be required to complete the satellite-to-satellite 
tracking surveys in these two regions. Again, some redundant information will 
be obtained, since the coverage regions around 94° W and 35° E overlap some- 
what. This will provide a valuable opportunity for cheeks and comparisons. 

Some satellite-to-satellite tracking will also be done while ATS-F is enroute 
between these 2 locations. This will, again, permit correlative studies which 
will strengthen the overall results. Perhaps another 50 hours or so could be 
devoted to this pahse of the activity. 

b. Orbit Determination Studies 

A satellite-to-satellite tracking and orbit determination investi- 
gation can involve data taking patterns such as those indicated in Figures 24 and 
25. 

These data sequences could presumably be abstracted from those indicated 
above for the scientific investigations. Accordingly, they would not necessarily 
involve the additional use of the GEOS-C /ATS-F satellite-to-satellite tracking 
resources. 
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REVOLUTION: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 


= ON E GEOS-C/ATS-F SATELLITE-TO-SATELLITE TRACKING 
PASS SPINNING HALF A REVOLUTION 

Figure 24. Satellite-To-Satellite Tracking Data for Orbit 
Determination Analyses 
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"o" S A ONE-MINUTE GEOS C/ATS-F SATELLITE-TO-SATELLITE 
TRACKING PASS. TRACKING PATTERNS SUCH AS THE FOLLOWING 
WOULD BE OF INTEREST: BD, BCD, ACE, & ABCDE. 

Figure 25. 
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B. Altimeter Data 


1. Introduction 

It was pointed out above that the altimeter will be useful for both 
gravimetric and oceanographic investigations and that, to a certain extent, these 
studies may be inter -twined. 

For example, the portion of the gravimetric geodesy investigation which is 
conducted with the aid of the altimeter will make use of the available knowledge 
of all of the factors which affect ocean height representation. These include, in 
addition to the gravitational field, the following oceanographic effects: 

1. Tides 

2. The General Circulation 

3. Currents 

4. Sea State 

5. Storm Surges 

6. Tsunamis 

The following discussion is concerned partly with the representation of these 
oceanographic effects using available information. It is anticipated that each 
effect represented will also be investigated with the aid of the altimeter data. 
Thus the above list of oceanographic topics is both a list of effects to be repre- 
sented, in connection with the gravimetric geodesy investigation, and a list of 
investigations to be conducted with the altimeter data. These various investi- 
gations are thus inter-related. In principle, with enough data, it should be 
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possible to separate the gravitational effects, which are nearly invariant over 
the lifetime of the experiment, from the oceanographic ones many of which vary 
with time in one way or another. In practice, however, it may be some time 
before enough data are available to make it feasible to effect such a separation 
without studying the data with both gravimetric and oceanographic thoughts in 
mind. In other words, it may be advisable that the gravimetric and oceanographic 
analyses be viewed as coordinated or joint undertakings in some sense. The 
observational as well as the analytical aspects of the gravimetric and oceano- 
graphic investigations may also be approached jointly. 

This portion of the discussion is aimed primarily at the gravimetric investi- 
gation and the representation of those oceanographic effects which are necessary 
in order to conduct the gravimetric investigation properly. Portions of the 
discussion that seemed to be especially likely to involve the possibility of joint 
planning such as those concerned with the observational programs, for example, 
do include suggestive thoughts relating to certain aspects of possible oceano- 
graphic investigations. 

2. Ocean Surface Altitude and Satellite Position Representation 

The representation of ocean surface altitudes and satellite positions 
will be considered in terms of the theoretical formulation, the organization of 
the calculations, and the specification of physical features associated with 
gravitational and oceanographic aspects of the earth. 
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a. Theoretical Formulation 


The theoretical formulation of the problem can be considered 
conveniently with the aid of the following definitions and relations : 
r (t) = spacecraft position vector 

h n (t) = vector normal to reference ellipsoid extending from reference ellipsoid to 
spacecraft 

h g = height of geoid above ellipsoid of reference 

h s = height of sea surface above geoid 

Bn(t) = r(t) - h n (t) 

h g (t) — h g h* (t) 

l^ (t) ^ h $ h* (t) 

h(t) = hj, (t) - h g (t) - h s (t) 

h (t) = I h (t) | 

Functional dependences can be indicated further as follows, where the 
symbols have the meanings specified. 

r(t) - r j t, e 0 , GM, C nm , S nm , 

P (h, <fi, X, t, ® a ,B a ), 

®, C, O, 7T, 7, 

p.m., UT1, T e | , 

where 


e 0 = (a, e, i, £2, w, m fl ), 


or 


-0 ~ Io’-0’ 
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GM GM = product of gravitational constant and mass of earth. 


C nm , S nm = cosine and sine coefficients of spherical harmonic potential terms. 

p (h, i ip, X, t» ® a , B a ) = atmospheric density as a function of: 

h: height 

if: latitude 

X: longitude 

t: time 

® a : solar activity 

B a : magnetic activity 
and where 

®: solar gravitational field 

c: lunar gravitational field 

a: radiation pressure 

tt: precession 

v : nutation 

p.m.: polar motion 

UT1 : earth’s rotational position 

T e : earth tides 

^g - hg I ^nm, ^nm’ ® ^ij» Nj)> 

**k’ g ^pg’ ^pg) | > 

where the 

^ij ) denote gravitational anomalies specified in terms of latitude and 
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longitude regions, the Mk denote gravitational anomalies specified in terms of 
mascon parameters, and the g(^ pq A pq ) denote gravitational anomalies specified 
either directly or in terms of equivalent geoid heights at a number of latitude and 
longitude points. 

Values at intermediate points are found by interpolation. A typical interval 
size would be specified as an integral multiple of a hundredth of a degree. The 
interval size may vary up to ten times over the range, and the interval patterns 
for latitude and longitude would be specified separately. Similar remarks apply 
to the specification of the order of interpolation. It is estimated that up to ten 
sets of the type g(sP pq A pq ) containing a total of 10 5 functional values might be 
needed initially, and that each of these storage parameters might be an order of 
magnitude larger for more advanced studies. 

The calculation of geoid height quantities, which are based on gravitational 
field expressions, could be implemented in more than one way. It could be done, 
for example, by evaluating these expressions at each point of interest. It could, 
alternatively, be done by evaluating these expressions at a specified set of times, 
e.g. , at the times corresponding to those employed in the numerical integration 
processes, and by employing numerical procedures to calculate the geoid height 
quantities at the points of interest. The method could be optimized keeping in 
mind numerical analysis and computational efficiency considerations. 
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One or more specified sets of the quantities g in (i£ pq) X pq )* e*g*» those 
associated with features such as trenches, say, would be multiplied by factors, 
k m , respectively, which would initially have the value unity, but would be solved 
for as unknown parameters. Ten to a hundred such quantities could be of interest. 
Thus, e.g. , 

n 

hg ~ ^6 | ^nm’ ® ^ij’ Np’ j + ®m ^pq’ ^pq^‘ 

m = 1 

_ | ^0’ ^0 ^pq» \>q^ | * 

where T 0 denotes ocean tides and the h 0 (y? pq , X pq ) denote heights above the 
geoid of oceanographic and meteorological features, such as currents, tsun am is, 
and wind and pressure fields, specified at a number of latitude and longitude 
points. The remarks about interpolation, intervals, arrays and unknowns, k m , 
made in connection with the functions g(ip pq , X pq ), etc. , apply in analogous 
fashion to the functions h 0 (<p pq , X pq ), etc. Other functions may also be used to 
represent oceanographic and meteorological phenomena. R n (t) = R n (t,R e ,f), 
where R e denotes the mean equatorial radius, and f denotes the flattening. 

The quantities r, h g , h s , and R^, will in general, be functions of unknown 
parameters which can be identified or associated with the indicated arguments. 
These unknown and any others involved, will be denoted here simply by x i . The 
partial differential coefficients can be found with the aid of methods and expres- 
sions employed in systems such as the Goddard Trajectory Determination System 
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(GTDS) , the Goddard General Orbit Determination System, Geostar, etc. That 
which is obtainable involves the use of numerical differentiation, numerical 

integration, or analytic theory, as appropriate. 

3h dr dRjj 3h g 3h s 

9xj 3xj ~ dXj “ 3xj 3Xj 

Here , for example , 

9h g 

3k^ “ 8m ^pq’ *W’ 

for m = 1, 2, . . . , n. 

3h/3xj can be derived, for example, numerically by differencing over a 
suitable time interval, e.g. , one which corresponds to the measurement interval. 
The determination of the time to be associated with the derived quantity can be 
made by an appropriate method. Alternatively, such a measure may be regarded 
directly as a height difference over a finite time interval and treated accordingly. 

b. The Organization of the Calculations 

The organization of the partial derivative calculations is indicated 
in Figure 26 and Table XVII. The notation of Figure 26 corresponds to that 
employed in reference 37 . 

c. The Specification of Physical Features 

Analysis of the altimeter data involves the representation of the 
gravitational field and of the types of oceanographic effects which were listed 
above in the introduction to this section. 
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Figure 26, 
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Table XVII 
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It will be assumed for the purposes of this portion of the discussion that the 
GEOS-C altimeter will be capable of instrumental resolution at the 2 meter level 
in a global mode and at the 0.5 meter level in a intensive or local mode. 

i. The Geoid 

The broad features of the geoid are given by satellite analysis. 
Typical representations are seen in Figure 4. 

Detailed studies of the geoid in certain regions have been conducted by a 
number of investigators including Von Arx, Talwani, Uotila, Rapp, Strange, 
Vincent, Berry, and Marsh.' ( 42, 47-51 ) An indication of the availability of results 

in several regions of interest off the East and Gulf Coasts of the United States 
is presented in Figure 27. An example of such a detailed geoid, and a satellite 
geoid for the same region, are seen in Figures 28 and 29. ^ 52, ***) 

ii. Tides 

The departure of the ocean surface from the geoid can be thought 
of in terms of tidal phenomena and the features of the general circulation, in- 
cluding current patterns and more vaiable features associated with currents such 
as the meanders of the Gulf Stream. 

Tidal phenomena have been studied by Hendershott, Munk and Zetler^ 53,54 ). 
Hendershott and Munk describe a model involving a number of amphidromic 
points which can be obtained by solving the Laplace Tidal Equations (LTE) in 
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the case where coastal values are given. (Cf. Figure 30) Cotidal lines for the 
Atlantic Ocean are given in more detail in Figure 31 which is due to Hansen. (SS) 



Figure 30. Cotidal and corange lines for the M 2 tide obtained by solving LTE 
with coastal values specified. Cotidal lines radiate from amphidromic points, 
corresponding to the progression of tidal crests in the sense indicated by the 
arrows around these points of vanishing tidal range (Section 5). High tide occurs 
along the cotidal lines labeled 0° just as the moon pases over Greenwich meri- 
dian. Successive cotidal lines delineate tidal crests at lunar hourly intervals. 
(For clarity, only selected cotidal lines are labeled, 30° corresponds to a dealy 
of one lunar hour.) Corange lines (5, 10, 20, 30, 75, 100, 125cm) connect 
locations of equal tidal amplitude (not double amplitude) . They surround amphi- 
dromic points (where range vanishes and phases very rapidly) and range max- 
ima (where range is longest and phases nearly constant). 
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Theoretical tides of Atlantic Ocean. Full lines: co-tidal lines referred to moon- 
transition through meridian of Grw.. dashed lines: co-range lines of the semi-diurnal tide 
M t in m (according to Hansen). 

Figure 31. 
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Hendershott's recent research has extended them to include effects of solid 
earth tides which he found to be significant (S6 >. 

iii. The General Circulation of the Oceans 

Problems associated with the general circulation have been dis- 
cussed by H. Stommel and K. Bryan. (Cf. references 57 and 58.) 

A model for the world-wide general circulation has been presented by 
Stommel. It appears in Figure 32. More detailed studies of particular regions 
such as the North Atlantic are being conducted by Bryan and Holland (cf. refer- 
ence 77). 


iv. Currents 

The general circulation embraces a number of major current sys- 
tems including the Gulf Stream, the Kuroshio Current near Japan, and the 
Antarctic circumpolar current. 

The Gulf Stream is of special interest here for several reasons including 
the fact that it meanders. 

a . The Gulf Stream Meanders 

Considerable effort has been devoted to the study of the meanders 
of the Gulf Stream. This research is discussed by Hansen (59) . The Gulf Stream 
proper can be thought of as a Kelvin wave. Departures from this model are 
also of interest here. 
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Figure 32. Dynamic Topography of the Sea Surface 
Relative to the 1000 -Decibar Surface 


In the study described by Hansen, the position of the Gulf Stream between 
Cape Hatteras and approximately 60° west longitude was observed at intervals 
ranging from a few days to a month. The 15° isotherm was chosen as an index 
of the main thermal front at a depth of 200 meters, and hence as an indication 
of the Gulf Stream position 59. 

Wavelike features having wavelengths of 200-400 km were observed. Hansen 
shows them in Figure 33. The evolution of the meanders in terms of phase pro- 
gression as indicated in Figure 34. It is seen there that the meanders progress 
at a rate of the order of a mean wavelength, taken to be 300 km, each couple of 
months. This corresponds to about 5km/day. 
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Figure 33 
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Figure 34. 


Altitude variations across much of the Gulf Stream due to the geostrophic 
effect are of the order of a meter. 

Some information about its position will be available from sources other 
than GEOS-C. 

The inclination of 65° also makes it possible to study the antarctic circum- 
polar current which is of interest for more than one reason. It is the only con- 
tinuous current which corresponds to a closed circuit or curve on the globe. It 
is also the strongest current, being equivalent in terms of mass flux to five or 
six Gulf Streams. The 65° orbit would also provide good geometry for viewing 
the Gulf Stream. This point is discussed in Section m, 4, C. 
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v. Sea State 


It is expected that the problem of interpreting the range measures 
in the presence of sea state effects will be attacked by an approach such as the 
one offered by Pierson and his colleagues* 60 *. 

vi. Storm Surges 

Storm surges, such as those associated with hurricanes, for 
example, can have amplitudes of up to a couple of meters and wavelengths of 
some fifty kilometers. Near coast lines the amplitudes have been known to in- 
crease to the order of seven or eight meters. 

Phenomena such as storm surges will be represented and/or observed when 
and where they occur, to the extent that this is practicable. 

vii. Tsunamis 

Tsunamis are thought to have amplitudes of nearly half a meter 
and wavelengths of perhaps three hundred kilometers. Detailed theoretical repre- 
sentation of a tsunami involves the representation of the ocean bottom topography 
and complex ray-tracing calculations to map out the advance of the expanding 
wave. The occurence of a tsunami in a region and at a time when it might be 
observed by GEOS-C would be a relatively rare event. It could, however, be 
modeled and studied in this way. 
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3. The Calibration of the Altimeter 


a. Short-Are Tracking of GEOS-C in the Caribbean Area 
The consideration of short-arc and long-arc tracking error budgets 
can begin with a look at the overall error problem. A typical error breakdown 
for the GEOS-C altimeter is indicated in Table XVIII (61 i Quantities associated 
with factors other than the orbit errors have an rms value of approximately 3 
meters. This leaves 4 meters or so which might be assigned to the calibration 
process if the 5 meter rms overall accuracy goal is to be met. Allowing 1 or 2 
meters for uncertainties associated with the geoid means that the uncertainties 
associated with the orbit determination process should contribute no more than 
about 3.5 meters. 


Table XVm 

GEOS-C Mission Altimeter Evaluation 


Satellite Altimeter System 
Measurement Error Source 

Error 

(m) 

Altimeter Instrumentation 

2 

Refraction 

0.2 

Reflection from Waves 

0.5 

Spacecraft Attitude 

2 

Root Sum Square 

2.9 

Calibration Error Allocation 

4.1 

Altimeter System 

5 

Evaluation Goal 
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If comparable accuracies are to be achieved over extensive areas, the 
accuracy within a short-arc calibration region such as the one in the Caribbean 
must be increased still further. This point will be discussed further below. 

A detailed analysis of short-arc tracking using lasers and cameras in the 
Caribbean area has been conducted by Berbert and Loveless (62) . Their analysis 
of a number of cases involving various combinations of lasers and cameras is 
summarized in Figure 35. Assumptions underlying these analyses are listed in 
Table XIX. The other angle measure accuracies of 100” listed there were those 
assumed for the laser angles used in the analysis indicated by the open circles 
in Figure 35. In all cases, in addition to the altimeter bias uncertainty, un- 
certainties in orbital, survey, and range measure parameters were also esti- 
mated. The triangle corresponds to a similar analysis of a three-laser-only 
case made for a much larger triangle based on stations at Antigua, Key West, 
and Panama. It resulted in a value of 4.1 meters, only slightly higher than 
that for the smaller triangle. As can be seen, a number of cases meet both the 
basic 4 meter requirement and the 3. 5 meter figure obtained by allowing a 
couple of meters for uncertainties associated with the geoid. 

Berbert and Loveless concluded that the 2-laser, 2-camera combination was 
probably the most cost effective in terms of the probabilities of obtaining 
reasonable amounts of data. 
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Table XIX 


GEOS-C Mission Altimeter Evaluation Analysis 


Assumptions 

A Priori 
Uncertainties 

Noise 

rms 

Recovered Quantities 



Range Measures 

2 m 

2 m 

Altimeter Height Measures 

100 m 

10 m 

Orbit - R&V Components 

1 km, 1 km/sec 


Station Positions — E, N, V Components 

30, 30, 1 m 


Other Quantities 



Camera Angle Measures 


1" 

I Other Angle Measures 


100" 


Tracking using lasers having ten centimeter accuracies can yield more 
accurate determinations of the orbit as has been pointed out by Vonbun. (63) Re- 
sults for 30 cm systems are indicated in Figure 36. This basic tracking system 
accuracy should give relative altitude accuracies of the order of half a meter. 
The limiting factor in this case would then appear to be the state of knowledge 
of the geoid which is thought to be on the order of a meter or two. 

b. Long-Arc Tracking of GEOS-C 

The surveys of the gravitational field over longer arcs will be greatly 
facilitated by the long-arc satellite -to -satellite tracking of GEOS-C which can be 
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Figure 36. GEOS-C Height Errors Using Laser Tracking 



conducted through ATS-F. The accuracy capability of this tracking approach is 
indicated in Figure 37. In the case looked at here, accuracies of some four 
meters or better persisted for almost three hours beyond the time interval 
shown in the Figure before the results deteriorated. Error analysis simulation 
results thus indicate that altitude accuracies in the three to four meter range 
can be achieved in this way. This is reasonably comparable to the current 
estimates of the accuracy of the world-wide geoid obtained from the satellite 
orbit analyses (18) . Satellite contributions to these have a spatial resolution of 
the order of 18°, however. Hence, altimeter and satellite -to -satellite tracking 
surveys even at the 6° resolution level will definitely provide new information. 
They will of course also provide the extremely valuable independent views which 
are so important. 

4. Ocean Surface Altitude Representation and Analysis Using 
Altimeter Data 

The set of GEOS-C data will itself make possible a better representa- 
tion of the effects indicated above. More than one kind of study holds promise. 
Tide can be analyzed on a global basis using all the altimeter data. Information 
about the geoid and the sea state can be sought everywhere. Current phenomena, 
such as those associated with the Gulf Stream meanders, will be watched in 
certain regions of the oceans. Storm surges and Tsunamis will be looked at 
when and where they occur, insofar as this is feasible. 
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Figure 37. 




a. Gravimetry 


Gravimetric investigations by means of a GEOS-C altimeter will 
include at least three kinds of studies. First, it is anticipated that it should be 
possible to conduct an altimeter survey having a spatial resolution of about 6 
degrees using the orbit ground track pattern which was worked out above in the 
discussion of the selection of a typical GEOS-C orbit. This will permit a direct 
inter-comparison with the results of the gravimetric investigation at the 6° res- 
olution level conducted by means of a satellite-to-satellite tracking system. A 
mutual validation of both the two approaches achieved in this way will increase 
the confidence of both types of results. 

A similar altimeter gravimetric survey could be conducted at a finer res- 
olution, e.g. , down to 1° . This can be thought of as the second type of experi- 
ment. This could also be conducted through the use of the GEOS-C orbit ground 
track pattern indicated earlier in the discussion of the GEOS-C orbit selection. 

Thirdly, fine structure could be investigated in greater detail in regions of 
special interest such as those near the East and Gulf coasts of America. 

The interesting region in the neighborhood of the Puerto Rican Trench, for 
example, could be studied with the aid of the Caribbean facilities indicated in 
the calibration discussion. 

Another possible region of tills type will be discussed in a later section. 
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The results of these regional studies, will, it is expected, be correlated 
closely with the surface gravimetry data discussed above. 

b. Tides 

Zetler has pointed out, for example, that GEOS-C will, over its life- 
time, provide data which can be used as the basis for a global tidal study. (54) 

As possible regional tidal study is discussed below in Section in, 4,C,i,B. 

c. A Region for Earth and Ocean Dynamics Studies 

As has been indicated above, surveys of certain phenonema such as 
the geoid and tides can be conducted on both global and regional bases. The 
following discussion deals with an area in which regional investigations of more 
than one type could be conducted fairly readily from a practical standpoint, 

i. Ocean Dynamics 

The Atlantic region off the coast of the Northeastern United States 
is of particular interest from the standpoint of the Gulf Stream meanders as was 
indicated in Figure 33 which is given by Hansen <59) . These features have ampli- 
tudes on the order of a meter. Tidal variations in this same region, while not 
quite as large as those found elsewhere, are nevertheless of considerable size, 
i.e. , of the order of a meter also. This is indicated in Figures 30 and 31 where 
certain tidal components are seen (53 & 55) . 
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This region is also a reasonably attractive one from the standpoint of some 
of the practicalities of short-arc tracking. Good advantage could be taken of 
lasers which are usually available at Goddard, Wallops Island, and possibly 
also at SAO. Other possible locations for lasers include Florida, Bermuda, the 
Canal Zone, and the Antilles, at Antigua and/or Grand Turk. C-band radars 
are also located at Wallops Island, Cape Kennedy, Bermuda, Antigua and Grand 
Turk. It would be valuable also to have an additional site near the point on the 
Bermuda meridian which is as far north of Goddard as Goddard is north of 
Bermuda. A laser at this location which is near the coast of the Bay of Fundy 
would be of interest not only from the point of view of oceanographic studies but 
also for earth dynamics studies, as will be pointed out later in the discussion. 

A site at this point, which is in Canada, or in nearby Maine would be suitable. 

It is anticipated that lasers will be located at some of the continental US sites 
shown in Figure 38, and that at least a couple of portable lasers will be available 
at times for use at some of the other sites indicated there although no such 
arrangements have been made at the present time. Three of these lasers in 
any of the triangular configurations in that map operating at the 10 cm accuracy 
level would make it possible to determine the altitude of GEOS-C with a relative 
accuracy of the order of half a meter or better in the neighborhood of the tri- 
angle. A fourth laser nearby would provide the important checks on the instru- 
mental biases by furnishing the redundant information. It would also reduce the 
impact of the cloud cover problem. 
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a. Gulf Stream Meander Studies 


The altimeter tracking patterns are also good for observing the 
Gulf Stream meanders in this region. Shown in Figure 39 are surface tracks of 
the 65° orbit with 6. 3° daily spacing which was obtained in the earlier discussion. 
It is seen that the northward and sourthward going tracks cross the two principal 
branches of a typical Gulf Stream meander nearly orthogonally, providing almost 
ideal geometry for studying the behavior of these interesting features. Each 
ground track seen in Figure 39 will be followed four days later by one removed 
just one degree from it, hence it will be possible to observe each feature once 
every four days. This frequency is well matched to the observational needs of 
a Gulf Stream meander experiment, as can be seen from inspection of Figures 
34 and 39 (59) . The mean wave length of ameander is often of the order of 300 kilom- 
eters, as was indicated above in the discussion associated with Figures 33 and 
34. A typical meander moves a distance equal to its own wave length in about 
a couple of months. This interval might be thought of as a characteristic time 
constant which can be associated with the Gulf Stream meanders in this sense. 
Observations every four days are well suited for a Gulf Stream meander experi- 
ment. In fact observations every ten days or so would be most welcome, as 
Hansen has already pointed out (59) . This also allows a margin for gaps in the 
observing program which might be due to such things as weather conditions or 
operational factors. 
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A CANADA 



Similar studies of the Kuroshio current could be conducted by means of 
lasers similarly placed in Japan and nearby islands such as IWO, Jima, say. 

B. Tidal Studies 

1. An Anti-Amphidromic Region Tidal Study 
A model for the M-, semi-diurnal component of the tide due to 
Hendershott appears in Figure 30. Hendershott has pointed out that the regions 
separating the amphidromic points in which the tidal range is large and the phase 
is slowly varying are of greater interest than the amphidromes themselves be- 
cause the tidal potential energy is largely to be found in these anti -amphidromic 
regions. He suggested that the most important of these regions have their cen- 
ters at approximately 0°N, 15° W; 15° S, 75°E;5°N, 135°W; 10°S, 115°W. (78) 
Zetler has also called attention to the importance of these regions (79) . 
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These locations are indicated in Figure 40. The first of these regions, 
which is the Sierra Leone Basin Area, is perhaps potentially the most suitable 
for study by means of the GEOS-C altimeter and orbit determination system since 
it is adjacent to island and coastal sites at which tracking lasers might be located, 
A US tracking facility is on Ascension Island. The French already have had a 
laser at Dakar and operated ground terminal facilities at Quagadougou. The 
§yncom satellites were tracked by a radio range and range rate system from a 
ship at anchor in the harbor of Lagos, Nigeria. These would be suitable sites 
from the scientific standpoint for tracking by means of 10 cm lasers. A set of 
three lasers at Freetown, Abidjan and Ascension might be even somewhat better 
from the geometric standpoint. No arrangements have been made at the present 
time for 10 cm laser tracking from any of these 5 sites. The location in the 
Lagos/Abidjan area will be referred to simply as the Gulf of Guinea site. 

Tracks of GEOS-C in this Sierra Leone Basin anti-amphidrome laser tri- 
angle region are seen in Figure 41. Twice each day the ground track passes 
through or very close to this triangle. At least one northbound track of the type 
seen, for example, would occur each day in the neighborhood of the triangle. 
These tracks are seen to be nearly orthogonal to the co-range lines of the semi- 
diurnal tide. Southbound tracks, nearly parallel to the co-range lines also occur 
daily in useful locations. The orbit worked out for GEOS-C in the above dis- 
cussion has the property of moving about 10.5 degrees each day relative to the 
moon. A complete cycle of the semi-diurnal lunar tide can be observed by 
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Fig. 40. COTIDAL AND CORANGE LINES FOR THE M z TIDE OBTAINED BY SOL- 
VING LTE WITH COASTAL VALUES SPECIFIED. COTIDAL LINES RADIATE FROM 
AMPHIDROMIC POINTS, CORRESPONDING TO THE PROGRESSION OF TIOAL 
CRESTS IN THE SENSE INDICATED BY THE HEAVY ARROWS AROUND THESE 
POINTS OF VANISHING TIDAL RANGE. HIGH TIDE OCCURS ALONG THE COTIDAL 
LINES LABELED 0° JUST AS THE MOON PASSES OVER GREENWICH MERIDIAN. 
SUCCESSIVE COTIDAL LINES DELINEATE TIDAL CRESTS AT LUNAR HOURLY IN- 
TERVALS. (FOR CLARITY, ONLY SELECTED COTIDAL LINES ARE LABELED,30° 
CORRESPONDS TO A DELAY OF ONE LUNAR HOUR.) CORANGE LINES (5, 10, 20, 
50, 75, 100, 125 cm) CONNECT LOCATIONS OF EQUAL TIDAL AMPLITUDE (NOT 
DOUBLE AMPLITUDE). THEY SURROUND AMPHIDROMIC POINTS (WHERE RANGE 
VANISHES AND PHASES VARY RAPIDLY) AND RANGE MAXIMA (WHERE RANGE IS 
LARGEST AND PHASES NEARLY CONSTANT). 

O SITES AT WHICH LASERS MIGHT BE LOCATED 
FOR LIMITED PERIODS 

• ATS-F SUBSATELLITE POINT 


□ APPROXIMATE CENTER OF ANTI-AMPHIDROM 1C REGION 
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Fig.41 T HEORETICAL TIDES OF ATLANTIC OCEAN. FULL 
LINES: CO-TIDAL LINES REFERRED TO MOON-TRANSITION 
THROUGH MERIDIAN OF Grw., DASHED LINES: CO-RANGE 
LINES OF THE SEMI-DIURNAL TIDE M 2 IN m (ACCORDING TO 
HANSEN). 

□ APPROXIMATE CENTER OF ANTI-AMPHIDROMIC REGION 

• SITES AT WHICH LASERS MIGHT BE LOCATED FOR 
LIMITED PERIODS 
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The ATS-F satellite will pass directly over the Sierra Leone Basin anti- 
amphidrome region as it moves from 94° west longitude to 34° east longitude. 

It will be within about 30° of the nominal center of the region on the equator at 
15° west longitude for the duration of the 17 -day interval in which GEOS-C can 
sense a complete cycle of the semidiurnal lunar tidal component. Satellite-to- 
satellite tracking between ATS-F and GEOS-C would certainly augment the laser 
capability. It may make it possible to conduct a useful investigation with a 
smaller number of lasers. Its value would be enhanced if it remained within 
about 13° of the center of the region for the entire 17 day cycle. Such a drift 
rate in this region, which is about half the one currently planned, would pro- 
vide an observing interval of some 26° . This is comparable with the charac- 
teristic GEOS-C observing span which is associated with the fact that its equator 
crossings in the neighborhood of a given point from day to day range over a 26° 
longitude interval. 

The GEOS-C altimeter operating in the 0. 5 meter accuracy mode with the 
relative orbital altitude determined to about half a meter should permit tidal 
observations to about a meter or better. The GEOS-C systems thus offer the 
possibility of conducting an initial, exploratory tidal experiment in the Sierra 
Leone Basin anti-amphidrome region. 


Satellite -to -satellite tracking from ATS-F at 94° W together with laser 
tracking from American and Pacific sites would open up the possibility of study- 
ing the anti-amphidrome regions in the Pacific as is indicated in Figure 42. 

/V/ 


r 
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Fig -^.THEORETICAL TIDES OF ATLANTIC OCEAN. FULL 
LINES: CO-TIDAL LINES REFERRED TO MOON TRANSITION 
THROUGH MERIDIAN OF Grw. r DASHED LINES: CO-RANGE 
LINES OF THE SEMI-DIURNAL TIDE M 2 IN m (ACCORDING TO 
HANSEN). 

• SITES AT WHICH LASERS MIGHT BE LOCATED FOR 
LIMITED PERIODS 

(§) ATS-F SUBSATELLITE POINT 

□ APPROXIMATE CENTERS OF ANTI-AMPHIDROMIC 
REGIONS 
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Fig. 42a. COTIDAL AND CORANGE LINES FOR THE M 2 TIDE OBTAINED BY SOL- 
VING LTE WITH COASTAL VALUES SPECIFIED. COTIDAL LINES RADIATE FROM 
AMPHIDROMIC POINTS, CORRESPONDING TO THE PROGRESSION OF TIDAL 
CRESTS IN THE SENSE INDICATED BY THE HEAVY ARROWS AROUND THESE 
POINTS OF VANISHING TIDAL RANGE. HIGH TIDE OCCURS ALONG THE COTIDAL 
LINES LABELED 0° JUST AS THE MOON PASSES OVER GREENWICH MERIDIAN. 
SUCCESSIVE COTIDAL LINES DELINEATE TIDAL CRESTS AT LUNAR HOURLY IN- 
TERVALS. (FOR CLARITY, ONLY SELECTED COTIDAL LINES ARE LABELED.30 0 
CORRESPONDS TO A DELAY OF ONE LUNAR HOUR.) CORANGE LINES (5, 10,20, 
50, 75, 100, 125 cm) CONNECT LOCATIONS OF EQUAL TIDAL AMPLITUDE (NOT 
DOUBLE AMPLITUDE). THEY SURROUND AMPHIDROMIC POINTS (WHERE RANGE 
VANISHES AND PHASES VARY RAPIDLY) AND RANGE MAXIMA (WHERE RANGE IS 
LARGEST AND PHASES NEARLY CONSTANT). 

O SITES AT WHICH LASERS MIGHT BE LOCATED 
FOR LIMITED PERIODS 

• ATS-F SUBSATELLITE POINT 


□ APPROXIMATE CENTER OF ANTI-AMPHIDROMIC REGION 
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2. A Continental Shelf Tidal Study 


A tidal experiment in the neighborhood of a continental shelf could be 
conducted in the western North Atlantic in the region covered by the tracking 
sites indicated in Figure 38. Amplitudes of the order of a meter are indicated, 
for example, in Figures 31 and 42. It is anticipated that lasers having 10 centi- 
meter accuracy will be located at sites in this region in connection with the 
altimeter calibration and the Gulf Stream studies outlined above. Data taken in 
connection with these activities would also be useful in connection with such a 
tidal study. Additional observations, specifically designed for the tidal investi- 
gation would also be needed. We see from Figure 42, that, again, the GEOS-C 
altimeter satellite ground track passes through or very close to the laser track- 
ing region once each day. At least one of the northbound tracks of the type seen 
in Figure 42 for example, would occur every day. These tracks are, again, 
nearly orthogonal to the co-range lines of the semi-diurnal tide as is also 
apparent from the figure. Southbound tracks, nearly parallel to the co-range 
lines, also occur daily in useful locations in similar fashion. The GEOS-C orbit 
will, again, sample a whole semi-diurnal cycle in about two and a half weeks. 

3 . A Bay of Fundy Tidal Study 

Tidal amplitudes in the Bay of Fundy region are the largest on earth. 
Ranges as high as fifteen meters are found. They should be amendable to direct 
measurement by the GEOS-C altimeter system. They should provide what may 
well be the first opportunity for an actual tidal observation by the GEOS-C. They 
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Fig. 42b. THEORETICAL TIDES OF ATLANTIC OCEAN. FULL 
LINES: CO-TIDAL LINES REFERRED TO MOON-TRANSITION 
THROUGH MERIDIAN OF Grw., DASHED LINES: CO- RANGE 
LINES OF THE SEMI-DIURNAL TIDE M 2 IN m (ACCORDING TO 
HANSEN). 

• SITES AT WHICH LASERS MIGHT BE LOCATED FOR 
LIMITED PERIODS 
@) ATS-F SUBSATELLITE POINT 
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4, Global Tidal Studies 


The use of altimeter data gathered over extended periods of time has 
been studied by Zetler (52) . He finds the prospects to be promising. 

ii. Earth Dynamics 

The laser tracking regions of Figures 38 and 41 also would have 
other uses in connection with the Earth Dynamics side of the Earth and Ocean 
Dynamic Satellite Applications program (4) . 

a. Gravimetric Fine Structure 

Fine structure in the gravity field should be deducible from 
observations made in the general area covered in Figure 38, but perhaps some- 
what away from the immediate neighborhood of the Gulf Stream meanders. 

b. Polar Motions and Earth’s Rotational Rate Variations 

The leg from Bermuda to the Bay of Fundy area would be 

suitable for observing polar motion in the manner of the experiment conducted 
by Smith ^ 64 \ The links from Goddard to such a pair of sites taken together would 
also be useful for a companion experiment to observe the variations of the 
earth's rotational rate. Similar remarks apply to other sets of lines in Figure 
38 which are roughly in latitude directions. 

The Sierra Leone Basin laser triangle would also be ideal for monitoring 
variations in the earth’s rotation rate since it is astride the equator. 
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5. Altimeter Data Requirements 

a. Gravitational Field Surveys 

The altimeter survey conducted at the 6° resolution level would 
involve, as indicated above, portions of the equivalent of some 4 days of orbit 
tracks. Roughly speaking, some three days of "sea tracks" will occur, since 
the altimeter will be over land about a quarter of the time. Again, from the 
standpoint of 6 degree coverage, redundant information will be obtained near 
the maximum latitudes. As before, however, continuous coverage will be use- 
ful, at least to begin with. Accordingly, about 75 hours of altimeter observing 
would be required for the 6° survey, and some 450 hours for the one degree 
survey. 

Such a set of observing programs would provide the basic data for gravi- 
metric geodesy studies leading to the shape of the geoid, for example. 

b. Tidal Analyses 

These data will also be of real value in connection with a variety of 
oceanographic investigations including global tidal studies. 

Additional types of observational data sets which would be useful for tidal 
studies can be considered in terms of the earlier discussion associated with 
Figures 41 and 42 say, and the thinking about the GEOS-C orbit. 
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As pointed out there, the moon moves some 10. 5 degrees per day relative 
to the GEOS-C orbit of the type contemplated, and hence a cycle of the semi- 
diurnal lunar tidal component could be observed in about 17 days. 

The durations of the northbound and southbound tracks of interest here 
would each be in the 10-minute range. Hence, about twenty minutes a day of 
observing, or a total of some six hours in all would give data for a tidal study 
of this type. Repetitions of this type of study in different regions such as the 
Pacific and Indian Oceans and at different times would be valuable. In other 
situations, as many as five passes per day might prove to be useful. In such 
cases a total of about fifteen hours would be required for a study. Some half 
dozen such tidal experiments would involve approximately a hundred hours of 
altimeter observing time. Much of this might be in addition to the 75 or so hours 
for the 6° geoid survey. It would be included in the 450 hours or so needed for 
the 1° geodetic altimeter survey. 

c. Gulf Stream Studies 

The altimeter observation tracks described above in connection with 
the regional tidal study in the Western Atlantic would also serve as the basis for 
a study of the Gulf Stream meanders. 

From the standpoint of the geoid study, the 6° and 1° resolution observations 
could be obtained anytime, in the sense that the gravity field is invariant over 
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the lifetime, and to within the accuracy, of the GEOS-C experiment. Studies of 
tides and the Gulf Stream meanders require observations made in specific se- 
quences as indicated in the above discussion. Hence, these should be given 
first priority in the scheduling except, of course, for observations of short- 
lived phenomena such as those associated with tsunamis and storm surges. Ob- 
servations designed to fill in the remaining regions required to complete the 
gravity surveys could then be scheduled so as to utilize the next increment of 
spacecraft capability. 

d. Calibration 

The altimeter calibration, which would precede studies discussed 
above, could also be conducted in the laser tracking region of Figure 38 in areas 
away from the meanders, near the coasts where ground truth is available. The 
calibration passes might in time be combined with the program to observe tides 
and meanders. Calibration and validation may thus involve only a small utili- 
zation of the GEOS-C resources relative to that which is envisioned for the 
scientific investigations, per se. 

Satellite-to-satellite tracking may also be useful when combined with pre- 
cision laser tracking in the laser tracking region of Figure 38 in making obser- 
vations in the neighborhood of the amphidromic point in the North Atlantic seen 
in Figures 30 and 31. Such a region could also be a good one in which to make 
the cross over point checks which have been proposed by Stanley (61) . 
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IV. A SET OF SATELLITE-TO-SATELLITE TRACKING STUDIES 

A. Introduction 

It is presently planned to conduct satellite-to-satellite tracking experiments 
with GEOS-C and ATS-F and with NIMBUS-E and ATS-F. Thought is also being 
given to the possibility of equipping Atmosphere Explorer and Small Astronomy 
Satellite (SAS-C) spacecraft with the capability for conducting satellite-to- 
satellite tracking operations with ATS-F. Such capabilities would be of value 
from the standpoint of the gravimetric geodesy studies, and, in the case of the 
Atmosphere Explorer spacecraft, for example, for operational reasons as well 
as is indicated in references 65 through 67. The present discussion will con- 
sider in more detail the aspects having to do with the gravitational field studies. 
These possibilities are only in the thinking stage at the moment and in fact, it 
appears as this material goes to press that, for practical reasons, it may not 
be feasible to equip the AE-C spacecraft with satellite-to-satellite tracking 
capability. The discussion of these cases involving AE and SAS spacecraft is 
included here, though, for the purpose of illustrating the possibilities inherent 
in orbits such as these. 

It is presently estimated that the ATS-GEOS and ATS-NIMBUS tracking 
systems will have accuracies of about 7 meters in range and 0.07 centimeters 
per second in range rate for a ten second integration interval. It will be assumed 
for the purposes of this discussion that the links between AE and ATS and SAS-C 
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and ATS would have similar characteristics. The AE and SAS-C systems may 
turn out to be slightly less accurate, perhaps by a factor of two, than the NIMBUS 
and GEOS systems for tracking through ATS. For this discussion however, it 
is convenient to assume that all these tracking systems will be similar, since 
this facilitates the understanding and comparison of the potential contributions 
of the different missions. 

Estimated orbital parameters of the low altitude satellites involved are 
given in Table XX. The values shown for the altitudes of Atmosphere Explorers 
when they are in their nearly-circular, low-altitude orbits are used for planning 
purposes. Decisions will probably be made during each flight mission concerning 
the different height ranges which will actually be surveyed from the circular 
orbits in the 700 to 250 kilometer altitude range toward the end of the space- 
craft's active lifetime. Nevertheless certainkinds of observations can be made 
about an ensemble of orbits of the types seen in Table XX. 

First, the array of different inclinations should be of real value. For ex- 
ample, range rate data will be taken over a given feature by satellites traversing 
at different angles in the orbits having different inclinations. It is anticipated 
that observable effects will frequently occur in two such distinct orbital paths 
when they pass over a given geographical region. Evidence of this type will be 
helpful in sorting out real physical effects from others which may be associated 
in some way with the measurement and analytical processes. This point is dis- 
cussed in reference 68, and in reference 44 in connection with Figure 20. It is 
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Table XX 


Orbit Parameters Planned or Considered for Spacecraft Proposed for 


Satellite-to-Satellite Tracking Experiments 


Spacecraft 

Altitude 

(kilometers) 

Inclination 

(degrees) 

NIMBUS-E 

1110, circular 

100 

GEOS-C 

1000, circular 

115 

SAS-C 

550, circular 

3 

AE-C 

700, circular 

65 


600, circular 

65 


500, circular 

65 


400, circular 

65 


300, circular 

65 


250, circular 

65 


150, elliptical* 

65 


120, elliptical* 

65 

AE-D 

700, circular 

100 


600, circular 

100 


500, circular 

100 


400, circular 

100 


300, circular 

100 


250, circular 

100 


150, elliptical* 

100 


120, elliptical* 

100 

AE-E 

700, circular 

20 


600, circular 

20 


500, circular 

20 


400, circular 

20 


300, circular 

20 


250, circular 

20 


150, elliptical* 

20 


120, elliptical* 

20 


*Perigee altitude. 

Apogee altitude ~4000km 
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seen that a variety of inclinations will be available ranging from within 3° of the 
equator to within 10° of the pole. 

When satellites pass over regions having the same ground track but at dif- 
ferent altitudes, it may be possible to learn more about the anomalous regions. 
For example, information may be obtained about the horizontal extent and/or 
the depth of the features involved. This point is discussed further in reference 
44 in connection with Figure 20. 

One case is of special interest here. NIMBUS-E and the second Atmosphere 
Explorer, AE-D, will both have the 100° inclination, but will orbit at widely 
different altitudes. 

The spatial resolutions obtainable with these spacecraft in the sense of the 
discussion of Section III, A, 1, associated with Figure 21 are indicated in Table 
XXI. Also shown there are the values for the resolution nodal longitude interval, 
X , i.e. , the spacing between equator crossings which corresponds to the spatial 
resolution, r, and the inclination, i, in the manner indicated in Section III, A, 2, 
i.e. , X = r esc i. For a satellite orbit at a 40° inclination, but otherwise having 
the same characteristics as those indicated for the GEOS-C call in the Table, 
the resolution nodal longitude interval would be about 8. 9° for example. 

Gravitational fields which have been derived recently on the basis of a com- 
bination of satellite data and gravimetry have spatial resolutions of 11° or 
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Table XXI 

Orbit and Gravimetric Geodesy Experiment Parameters for Spacecraft 


Proposed for Satellite -to -Satellite Tracking 


Spacecraft 

Altitude 

(kilometers) 

Inclination 

(degrees) 

Spatial 

Resolution 

(degrees) 

Resolution 

Nodal 

Longitude 

Interval 

(degrees) 

NIMBUS-E 

1110 

100 

6.2 

6.3 

GEOS-C 

1000 

115 

5.7 

6.3 

SAS-C 

550 

3 

4.0 

90.0 

AE-C 

700 

65 

4.5 

5.0 


600 

65 

4.0 

4.5 


500 

65 

3.5 

4.0 


400 

65 

3.0 

3.4 


300 

65 

2.5 

2.7 


250 

65 

2.0 

2.4 


150* 

65 

1.5 

1.6 


120* 

65 

1.2 

1.4 

AE-D 

700 

100 

4.5 

4.6 


600 

100 

4.0 

4.2 


500 

100 

3.5 

3.6 


400 

100 

3.0 

3.1 


300 

100 

2.5 

2.6 


250 

100 

2.0 

2.2 


150* 

100 

1.5 

1.5 


120* 

100 

1.2 

1.3 

AE-E 

700 

20 

4.5 

13.2 


600 

20 

4.0 

12.0 


500 

20 

3.5 

10.2 


400 

20 

3.0 

9.0 


300 

20 

2.5 

7.3 


250 

20 

2.0 

6.3 


150* 

20 

1.5 

4.3 


120* 

20 

1.2 

3.6 


♦Perigee Altitude 
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so, or about 1200 kilometers near the equator. The geoids associated with them 
are considered to be reliable to about three meters. It is estimated that the 
information content of these fields which corresponds to spatial resolutions finer 
than about 18° , or about 2000 kilometers near the equator, is derived largely 
from the surface gravity data. <Cf. reference 18.) 

It is seen from Table XXI that the satellite -to -satellite tracking experiments 
which are contemplated offer the prospect of improving the spatial resolution by 
factors in the range from about half to one order of magnitude. 

The geoids associated with fields such as 1 through 5 in Table I have been 
estimated to have accuracies on the order of 15 meters (69) . 

The SAO 69 (II) solution in Table I, when truncated as 8, 8 does not repre- 
sent an improvement over its predecessor, the SAO M-l field. Except for 
resonant terms, the improvement due to the added satellite data in the latest 
solution is thus probably to be found chiefly in the terms of the 9th and 10th 
degrees, some of which are poorly determined (18) . The geoid associated with 
the general portion of that solution, as it reflects satellite data, is estimated to 
be somewhat better than fifteen meters but perhaps not greatly so. 

o 

Comparisons of this field with results obtained using recent gravimetric 
data indicate that geoid uncertainties of five to seven meters are to be expected 
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when the gravimetric data are moderately dense and accurate. (70) This field as 
it reflects satellite data per se, is probably still less accurate. 

It is estimated, thus, that the contribution of the satellite data to this field 
corresponds to a geoid accuracy on the order of ten meters, say, which in turn 
corresponds to some two and a halt milligals. Accuracies comparable to this 
should be obtainable from 700 kilometer altitude orbits with the ATS satellite- 
to-satellite tracking system using 45-second integration times to achieve accu- 
racies of the order of 0.15 to 0.2 millimeters per second. (Cf. Figure 22 and 
Reference 45.) A one-minute integration time would give a tracking capability 
on the order of 0.1 millimeters per second which corresponds to an acceler- 
ation resolution capability of the order of 1.5 milligals in terms of sensing mean 
anomalies in squares about 5° on a side for a satellite at an altitude of 700 kilom- 
eters. (45) In a sense, the ultimate range rate integration interval would corres- 
pond to the time required to pass over or traverse the square. A further im- 
provement in tracking accuracy and the corresponding sensing capability might 
be achieved in this way. On this basis, the Atmosphere Explorer, in an orbit • 
at about 250 kilometers altitude, would have an acceleration resolution capability 
of the order of a milligal. This may be somewhat optimistic, hence it will be 
assumed, more conservatively, that the acceleration resolution capability of the 
Atmosphere Explorer in a 250 kilometer orbit will be on the order of a couple of 
milligals. The acceleration resolution capability of GEOS-C and NIMBUS-E is 
estimated to be on the order of four milligals, and the acceleration resolution 
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of SAS-C is estimated to be on the order of three milligals. In making these 
various estimates, it was assumed in accordance with Schwarz's finding that 
the sensitivities vary roughly as the reciprocal of the altitude in this general 
range. (45> 71) 

The estimates of the acceleration resolutions to be derivable from satellite- 
to-satellite tracking appear to be comparable with the estimates of the accuracies 
we have now from satellites, however, the differences between these two kinds 
of estimates may be no more than their uncertainties. 

Certain additional features of interest which are associated with these dif- 
ferent contemplated experiments are indicated in the following discussion. 

B. NIMBUS-E 

The NIMBUS-E/ATS-F experiment is expected to be the first of the satellite- 
to-satellite tracking experiments. As was indicated in the above discussion, it 
will provide us with an improvement of nearly a factor of three in spatial reso- 
lution. It will also, of course, give us our first experience with this new type 
of data and technique. 

The 6.6 bidaily nodal longitude interval which is indicated in Table XXII 
implies that some 27 orbital arcs of interest will be trackable from ATS-F in 
a given position. There is both a northbound and a southbound pass associated 
with each interval hence some 54 arcs will be of interest. The total tracking 
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interval would thus be approximately 54 hours. ATS will at first be at 94° west 
longitude and later at 35° east longitude. A second survey centered at the 
eastern longitude would involve another 54 arcs and another fifty -four hours. 
Some fifteen hours would involve overlapping coverage which should be useful 
for correlative and corroborative purposes. In addition, as ATS moves from 
one location to another, satellite-to-satellite tracking could provide additional 
data having different geometrical characteristics. In all some three such sur- 
veys, centered at about 94° west, 35° east, and 30° west, would be very 
valuable. 


Table XXn 

Gravimetric Geodesy Experiment Parameters and Observation 
Requirements for Certain Spacecraft Proposed for 
Satellite-to-Satellite Tracking 


Spacecraft 

Altitude 

(kilometers) 

Resolution 

Nodal 

Longitude 

Interval 

(degrees) 

Daily or 
Bidaily 
Nodal 
Longitude 
Interval 
(degrees) 

Number of 
Arcs 

Per Survey 

Number of 
Hours Per 
Survey 

NIMBUS-E 

1110 

6.3 

6.6* 

54 

54 

GEOS-C 

1000 

6.3 

6.3 

57 

56 

SAS-C 

550 

90.0 

— 

4 

3 


*Bidaily 


106 







C. GEOS-C 

The GEOS-C satellite orbit is presently envisioned to be somewhat lower 
than the NIMBUS-E orbit. This will offer a corresponding increase in resolution. 
Also the GEOS-C orbit will not be subject to perturbations by control jets. In 
addition, the orbit of GEOS-C will be known much more accurately, independ- 
ently of the satellite -to -satellite tracking, through the use of the very accurate 
geodetic tracking systems which it will employ. Hence, it is the ideal satellite 
on which to really evaluate the satellite-to- satellite tracking system as was 
pointed out in reference 72 in which this experiment was first proposed. The 
GEOS-C satellite will, in addition, as was pointed out above, also provide data 
at another inclination. Observational requirements worked out along the lines 
indicated in the discussion of the NIMBUS-E case are indicated in Table XXII. 

D. SAS-C 

The SAS-C spacecraft, orbiting at a relatively low altitude, would afford a 
still further increase in spatial resolution and a gain of nearly a factor of two 
in acceleration resolution. SAS-C would also add data at still another inclination. 
Again, typical observational requirements are listed in Table XXII. 

E. The Atmosphere Explorers 

The Atmosphere Explorers would provide a marked increase in the capa- 
bilities for both spatial resolution and acceleration resolution. They would also, 
as was pointed above, provide data at different inclinations, and at a number of 
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heights for each of the inclinations. This would enhance considerably the value 
of such a gravimetric geodesy investigation. 

Gravimetric geodesy investigations which could be performed with Atmos- 
phere Explorer by means of satellite -to -satellite tracking are of two types, 
those associated with the elliptic orbits and those associated with the circular 
orbits. Present tentative plans call for an elliptic orbit for AE-C having an 
inclination of 65° , an apogee height of about 4000 km, and a perigee height 
which will usually be about 150 km, and which will be lowered to about 120 km 
periodically. The circular orbits will be at several heights in the range from 
about 250 to 700 km. 

For the elliptic case the orbit having a perigee of about 150 Ion would be in 
existence long enough to permit a survey to be made. The altitude of about 150 
km corresponds to a spatial resolution of 1. 5 degrees in the sense of the dis- 
cussion of Section III, A, 1, which is associated with Figure 21. The planned 
orbit is close to one which will have a daily nodal separation of 1.6 degrees. 
This would be the appropriate spacing for the case in which the perigee is near 
the equator, which is one of the possibilities now under consideration. If the 
perigee is far from the equator, the survey could be completed with a corres- 
pondingly larger nodal spacing and with a correspondingly smaller number of 
passes than is indicated here. 
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A slight increase in the period of the orbit, of the order of a minute, and 
a corresponding increase in apogee height, of the order of a 100 km, would give 
a daily nodal spacing of 1,6 degrees. This orbit, then, would permit the making 
of a survey which would provide the coverage to correspond to the lowest altitude 
reached by the AE-C spacecraft in such a case. The altitude region between 
150 and 1110 km would be of interest. This would permit correlation of results 
with those obtained from all the satellites to be tracked from ATS, i.e. , 
NIMBUS-E, GEOS-C, and possibly SAS-C . The Atmosphere Explorer will be 
below 1110 km over a true anomaly range of about 142° , or for about 34 minutes 
per revolution. Some 100 arcs under ATS when it is at a given location •would 
suffice to complete a survey at the 1. 5 degree resolution level. Roughly speaking, 
one pass every other day during the eight months or so at AE-C is in the ellip- 
tical orbit would be sufficient. An additional eighty passes would complete the 
survey if ATS would move between 94° west and 35° east during this period. 

Thus a total of some 180 passes, or less than one pass per day on the average, 
would suffice for the entire longitude history of ATS-F. The corresponding 
tracking time requirements are about 57 hours for the 100 passes and 105 hours 
for the 180 passes. 

If the perigee is at the equator, each arc which is of interest could be ob- 
served from ATS during one continuous interval. If the perigee is at the pole, 
the arcs of interest would be observed in two equal portions, when the ATS is 
on either side of the AE-C orbit. When perigee occurs in an intermediate 
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position the arcs of interest would be observed from ATS in two unequal time 
intervals occurring when ATS is on either side of the AE orbit. The total ob- 
serving time would be about the same for all these cases, i.e. , about 105 hours. 

The dipping down of the perigee to 120 km altitude provides a spatial reso- 
lution that is even finer, i.e. , about 1.2°. In view of the limited time during 
which the perigee will be at these low altitudes, however, a complete survey 
would probably not be practical. Nevertheless the observation of portions of the 
orbit at this height from ATS at different longitudes would be of great interest. 

The circular orbits planned for the AE missions would also be of great 
value from the gravimetric geodesy standpoint. There is interest in matching 
the nodal spacing at the equator to the spatial resolution capability associated 
with the altitude. It is seen from Table XXIII that the desired resolution nodal 
longitude intervals do not correspond too well with the daily or bidaily nodal 
longitude intervals associated with the orbits at the nominal altitudes at the 
100 kilometer intervals. 

The nodal spacing could be matched to the desired spacings more closely 
by modifying somewhat the periods and altitudes listed in Table XXIH. The 
orbital altitude parameters in or near the range from 250 to 700 km which are 
close to those considered for the Atmosphere Explorer C and which would also 
yield the nodal separations that could form the basis for good gravimetric 
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geodesy experiments are indicated in Table XXIV. It is seen that only relatively 
small changes from the nominal orbits would suffice in most cases. 


Table XXHI 

Gravimetric Geodesy Experiment Parameters for Certain 
Possible Atmosphere Explorer-C Orbits 


Altitude (kilometers) 

Resolution Nodal 
Longitude 
Interval (degrees) 

Daily or Bidaily 
Nodal Longitude 
Interval 
(degrees) 


5.0 

1.2* 


4.5 

6.9 


4.0 

0.8 


3.4 

2.9* 


2.7 

6.9 

250 

2.4 

2.9 


*Bidaily 


The approximate values shown in Table XXIV are suggested for illustrative 
purposes. It is anticipated that further study of the matter of orbit selection 
from the standpoint of the atmospheric research requirements and the gravi- 
metric geodesy criteria could reveal orbit parameters which would meet the 
atmospheric experimenters needs and at the same time make it possible to 
conduct valuable gravimetric geodesy experiments. 
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Table XXIV 


Experiment Parameters and Observation Requirements for Certain 
Possible Atmosphere Explorer-C Orbits Which Would Be Valuable 
for the Gravimetric Geodesy Investigation 


Altitude 

(kilometers) 

Resolution 

Nodal 

Longitude 

Interval 

(degrees) 

Daily or 
Bidaily 
Nodal - 
Longitude 
Interval 
(degrees) 

Number of 
Arcs per 
Survey 

Number of 
Hours per 
Survey 

Number of 
Arcs per 
45° Survey 

Number of 
Hours per 
45° Survey 

725 

5.1 

5.1* 

63 

45 

35 

17 

570 

4.4 

4.4 

73 

50 

41 

19 

460 

3.7 

3.7 

87 

58 

49 

22 

405 

3.4 

3.4* 

94 

62 

53 

23 

240 

2.3 

2.3 

139 

88 

78 

33 

160- 

4100 

1.6 

1.6 

100 

57 

i 




* Bidaily 


It is seen from Table XXIV that on the order of 60 to 140 passes of nearly 
half an orbit in length would be of interest at the heights shown there. The AE 
satellite will probably spend only some 10 or 30 days at each of these altitudes, 
however. Thus, from two to fourteen passes per day would be required, which 
may be beyond the resources available. In such a case, a more limited survey 
could be attempted. 
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A useful survey could be conducted over the region lying within about 45° 
of the sub-ATS region. This would involve only a much smaller number of arcs, 
and each one would not be nearly so long as in the case of the complete survey. 
Hence only a fraction of the observing time would be required for such a limited 
survey. Such a survey would correspond reasonably well to the type of experi- 
ment which was conducted in connection with the discovery of the lunar mascons. 

The satellite -to -satellite tracking system of the type which might be used 
between AE and ATS will have an accuracy of about 0.07 cm per second for a 
ten second integration interval. This corresponds to an acceleration of about 
seven milligals. It is planned that the accelerometer to be used in connection 
with the atmospheric density experiments will operate over two resolution 
ranges while the experiment is in progress. These extend, respectively, from 
5 x 10" 4 to 10 ' 6 g, and from 10 -5 to 2 x 10~ 8 g. It appears, then, that the accel- 
erometer resolution would be adequate for the needs of such a gravimetric 
geodesy investigation. 

Some of the parameters of geodetic interest associated with the spacecraft 
discussed here are summarized in Table XXV. 

V. COMPANION EARTH AND OCEAN PHYSICS STUDIES 

Precise laser tracking of the type contemplated in connection with the dis- 
cussion associated with Figures 38 and 41 has also yielded interesting new 


113 



Table XXV 


Spacecraft 

Inclination 

(degrees) 

Altitude 

(kilometers) 

Approximate 

Spatial 

Resolution 

(degrees) 

Order of 
Acceleration 
Resolution 
(milligals) 

Existing set 
of Spacecraft 



18.0 

2.5 

NIMBUS-E 

100 

1111 

6.2 

2.0 

GEOS-C 

115 

1000 

5.7 

2.0 

AE-C 

65 

725 

4.6 

1.5 

SAS-C 

3 

550 

4.0 

1.5 

AE-C 

65 

240 

2.0 

1.0 


results which have exciting implications for the emerging Earth and Ocean 
Physics Applications Program (EOPAP) now in the planning stages at NASA. (4) 

A. Polar Motion 

The Preliminary Polar Motion Experiment, has, for example, yielded one 
component of the motion of the pole which agrees with the BIH smoothed curve 
to within a meter over a five-month period. (73) These results appear in Figure 
43. The basic data were obtained by tracking the BE-C spacecraft with an 
accurate laser located at the Goddard Space Flight Center. The significance of 
this result lies in the fact that the laser values for the latitude were obtained 
with six-hour data spans from a single observing station, whereas the BIH curve 
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VARIATION IN LATITUDE, SEC. OF ARC 



Figure 43. Variation in Latitude of Goddard Laser 


METERS 



is derived from the smoothing of five day averages of data obtained from some 
40 stations around the world. 

The Preliminary Polar Motion Experiment configuration actually involved 
a pair of lasers, the one at Goddard and a second one located on the same 
meridian about 400 kilometers to the north, at Seneca Lake, New York. 

B. Baseline Determination 
1. A Laser Measurement 

The data from these two lasers were used to measure the distance 
between them. Results of several determinations of this intersite baseline in- 
volving data spans ranging from about a week to a couple of months are indicated 
in Table XXVI. It is seen that a repeatability about 30 centimeters was obtained. 


Table XXVI 

Godlas-Senlas Chord Distance 


Solution 


Survey Difference 

1 

408699. 20 m 

+ 43 cm 

2 

408698.87 m 

+ 10 cm 

3 

408699.33 m 

+ 56 cm 

4 

408699.44 m 

+ 67 cm 

5 

408698.91 m 

+ 14 cm 

Ground Survey 

408698.77 m 
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The same observational material was used in still another way to determine 
a value for Love's tidal number, k 2 . 

2. VLBI Measurements 

Very Long Baseline Interferometry (VLBI) has also yielded baseline 
determination results with accuracies of the order of a meter or two. (74, 7S> 
More recent results indicate that higher precision can be expected in the meas- 
urement of baselines by means of VLBI. (75) 

C. The San Andreas Fault Experiment (SAFE) 

The baseline determination results indicated in Table XXVI were obtained 
with lasers having a 50 centimeter noise figure. It is planned to utilize better 
lasers, having accuracies of the order of 5 to 10 centimeters, to study the 
motion of the tectonic plates at the San Andreas Fault. It is hoped that this San 
Andreas Fault Experiment (SAFE) will result in the determination of the base- 
line length with a precision of 10 centimeters. * 76) 

VI. THE GEOPAUSE SATELLITE SYSTEM CONCEPT 

The GEOPAUSE satellite system concept offers promising possibilities for 
conducting Earth and Ocean Physics Applications investigations of several types 
at and below the decimeter level. These include a set of earthquake-related 
studies of fault motions and the earth's tidal, polar and rotational motions, as 
well as studies of the gravity field and the sea surface topography which should 
furnish basic information about mass and heat flow in the oceans. 
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The state of the orbit analysis art is presently at about the 10 meter level, 
or about two orders of magnitude away from the 10 cm range accuracy capability 
expected in the next couple of years or so. The realization of a 10 cm orbit 
analysis capability awaits the solution of four kinds of problems, namely, those 
involving orbit determination and the lack of sufficient knowledge of tracking 
system biases, the gravity field, and tracking system biases, the gravity field, 
and tracking station locations. 

The GEOPAUSE satellite system concept provides promising approaches in 
connection with all of these areas. A typical GEOPAUSE statellite orbit has a 
14 hour period, a mean height of about 4. 6 earth radii, and is nearly circular, 
polar, and normal to the ecliptic. At this height only a relatively few gravity 
terms have uncertainties corresponding to orbital perturbations above the decim- 
eter level. The orbit is, in this sense, at the geopotential boundary, i.e., the 
geopause. The few remaining environmental quantities which may be significant 
can be determined by means of orbit analyses and accelerometers. The 
GEOPAUSE satellite system also provides the tracking geometry and coverage 
needed for determining the orbit, the tracking system biases and the station 
locations. (Cf. Figures 44 and 45.) Studies indicate that the GEOPAUSE satel- 
lite, tracked with a 2 cm ranging system from nine NASA affiliated sites, can 
yield decimeter station location accuracies. Five or more fundamental stations 
well distributed in longitude can view GEOPAUSE over the North Pole as is 
indicated in Figure 46. This means not only that redundant data are available 
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for determining tracking system biases, but also that both components of the 
polar motion can be observed frequently. When tracking GEOPAUSE, the NASA 
Goddard network sites become a two-hemisphere configuration which is ideal 
for a number of earth physics applications such as the observation of the polar 
motion with a time resolution of a fraction of a day. 

GEOPAUSE also provides the basic capability for satellite-to-satellite 
tracking of drag-free satellites for mapping the gravity field and altimeter 
satellites for surveying the sea surface topography. Geometries are indicated 
in Figure 46. It has been estimated that satellite-to-satellite tracking at the 
0.03mm per second level can yield 10 cm accuracy in the determination of the 
geoid. (3) GEOPAUSE tracking a coplanar, drag-free satellite for two months to 
0,03mm per second accuracy should yield the geoid over the entire earth to 
decimeter accuracy with 2. 5° spatial resolution. Two GEOPAUSE satellites 
tracking a coplanar altimeter satellite can then yield ocean surface heights 
above the geoid with 7° spatial resolution every two weeks. These data will 
furnish basic boundary condition information about mass and heat flows in the 
oceans which are important in shaping weather and climate. 

The GEOPAUSE satellite system thus offers the promise of making key 
contributions to the study of the major areas of earth and ocean dynamics. 
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